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Abstract We review theoretical developments in studies of dense matter and its phase structure 
of relevance to compact stars. Observational data on compact stars, which can constrain the 
' properties of dense matter, are presented critically and interpreted. 
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1 INTRODUCTION 

Since their discovery almost 40 years ago, neutron stars have been recognized 
as promising laboratories for studying matter under extreme conditions. All the 
ambient conditions that characterize these objects tend to the extreme. The 
typical density inside a neutron star is comparable to the density inside nuclei, 
Pnuclear = 2.5 x 10 14 g cm -3 , which corresponds to a number density of baryons 
tib — 0.15 fm~ 3 . Under these conditions the nucleons are very degenerate, with 
a typical Fermi momentum kp ~ {2>-K 2 nB) 1 ^ in the range of 300-600 MeV. In the 
interior, the interaction energy of the nucleons is several times larger than their 
Fermi degeneracy energy Ep = kp/2m ~ 60-150 MeV. The baryon chemical 
potential, which is defined as the energy needed to introduce one unit of baryon 
number, can be 500-1000 MeV larger, at the center of the star, than the rest 
mass of the nucleon. There the internucleon distance becomes comparable to 
their intrinsic size of 1 fm. Under such extreme conditions several new states of 
matter, containing new high-energy degrees of freedom, may be favored. Novel 
hadronic phases - which contain, in addition to nucleons, a Bose condensate of 
pions or kaons, and/or hyperons - have been proposed and studied in some detail. 
There is also the exciting possibility that the core may contain deconfined quark 
matter, much like the quark-gluon plasma being probed in terrestrial experiments 
such as the Relativistic Heavy Ion Collider (RHIC). Although model descriptions 
of these various phases are still in rudimentary form, in several cases qualitative 
differences between the exotic phases and nuclear matter have been identified. 
A quantitative understanding of these trends and, in particular, how they could 
affect observable aspects of compact star structure and evolution is an exciting 
area of research. A primary purpose of this review is to convey some of this 
excitement in a pedagogic style to make it accessible to graduate students. This 
review is not intended to be a comprehensive account of all recent developments. 
Instead, we have chosen topics that serve (in our view) as pedagogic tools that 
help illustrate key concepts of relevance to dense matter in compact stars. 

Figure ^ illustrates a theoretician's view of a neutron star. The environment 
just outside the star, characterized by large magnetic (and electric) fields and 
temperatures, is extreme from plasma physics point of view. At the surface of the 
star we expect a very thin atmosphere composed of hydrogen and in some cases 
perhaps a mix of heavy elements, or even a condensed magnetic surface. This 
surface is of utmost observational importance because it allows us to determine 
the temperature, and possibly the radius, of the neutron star, an issue we describe 
in Section \l. 21 An envelope exists just below the atmosphere where matter is not 
yet fully degenerate and, with a thickness of a few tens of meters, it acts as a 
thermal insulator between the hot interior and the surface. The outer 500-1000 
meters of the star, its crust, contains nuclei, forming a lattice immersed in a 
quantum liquid of neutrons, most probably in a superfluid state (Figure ^ insert 
A). Owing to the rotation of the star, the neutron superfluid forms vortices, which 
very likely pin on the nuclei and participate in the glitches. As we go deeper 
into the crust, when approaching nuclear density there is a first-order phase 
transition from the inhomogeneous regime of the crust to the homogeneous core. 




Figure 1: Inside a Neutron Star. 



In this regime, termed the nuclear pasta (Figure ^ insert B), as density grows, 
nuclei are increasingly elongated (the spaghetti phase), then form bidimensional 
structures (lasagnas) - with the space between them still filled by the superfluid 
neutrons - untill the geometry inverts in the Swiss cheese phase, where bubbles of 
neutrons are immersed in the almost homogeneous neutron+proton liquid. The 
physical reasons for the peculiar aspect of this phase transition are elucidated in 
Section EH At a density of approximately 0.6 /Onuciean we enter the core where 
the neutron superfluid coexists with the proton fluid, which is probably a Type 
II superconductor that reacts to the neutron star magnetic field by confining it 
into fluxoids (Figure ^ insert C). A few kilometers below may be the inner core, 
marked in Figure ^ by "?". This question mark is the focus of this review. 

In Section |2] we present a selected set of observational constraints. These in- 
clude mass (Section 12.1(1 and radii (Section I2.2|) measurements that probe the 
equation of state (EoS), i.e., the high-energy properties of dense matter. Mea- 
surements of temperatures and thermal luminosities of isolated neutron stars and 
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accreting ones in binary systems are presented in Section r2.2l and Section \'2.'6l re- 
spectively. These observations provide information about the structure of the 
stars, i.e., the EoS, but also about the low-energy response functions, such as 
specific heat, transport properties, and neutrino emission, which probe the low- 
energy (keV-MeV) excitation spectrum. In Section |3] we discuss the essential 
physics of dense nuclear matter and describe the current state of model calcu- 
lations of the nuclear EoS and low-energy response. In Section 21 we describe 
theoretical models for phase transitions to novel states of matter such as Bose- 
condensed hadronic matter, hyperonic matter, and quark matter at supranuclear 
density. Section confronts theoretical models of neutron star structure and of 
evolution derived from the studies described in the two previous sections with the 
observational baggage acquired in Section [3 Finally, we offer some conclusions 
in Sectional 

2 OBSERVATIONAL CONSTRAINTS 

Among the more than 100 billion stars in the Milky Way, almost 1% of them 
(~ 10 8 -10 9 ) are expected to be compact stars. Nevertheless, the number of 
known compact stars is of the order of only a couple thousand, the majority 
being radio pulsars (more than 1500) and the rest comprising approximately 
300 accreting stars in binary systems, plus a few dozen single stars that do not 
show any pulsed radio emission and were discovered through their high-energy 
emission. 

2.1 Mass Measurements 

The masses of almost 40 neutron stars in binary systems have been measured, 
and the results are summarized in Figure [21 The major source of uncertainty 
usually is the inclination of the orbit with respect to the line of sight, and it 
is only in highly relativistic systems that general relativistic effects allow for a 
complete characterization of the orbit (^Q) . 

The binary systems exhibited in Figure |2] are divided into four classes. The 
binary pulsar systems contain a radio pulsar whose companion is either a neutron 
star (PSR+NS) or a white dwarf (PSR+WD). In the PSR J0737-3039A and B 
system, both neutron stars are detected as radio pulsars, whereas in the three 
systems with the companion marked by "?" a white dwarf could be present 
instead of a second neutron star. The other two classes are X-ray binaries in 
which the compact object is accreting matter from its companion. In a high-mass 
X-ray binary (HMXB) the companion is a massive star, M c > IOMq, whereas in 
a low-mass X-ray binary (LMXB) the companion's mass is below 1M@. 

In a HMXB, once the companion ends its life in a core collapse supernova, 
producing a second neutron star, the system may be disrupted or remain bound, 
in which case it will most probably emerge as a PSR+NS binary system. A 
LMXB system will end its life once accretion stops, the companion's mass having 
been severely reduced, and the system will emerge as a PSR+WD system. Given 
the short lifetime of massive stars, accretion in a HMXB lasts at most a million 
years, whereas in a LMXB it can last for billions of years. Because accretion is 
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Figure 2: Mass measurements of neutron stars in binary systems. Values for 
pulsar + neutron star (PSR+NS) and pulsar + white dwarf (PSR+WD) systems 
are from References and with updated measurements of PSR J0751+1807 
from Reference ©, PSR J1713+0747 from Reference ©, PSR J1748-2446I & J 
from Reference and PSR J1909-3744 from Reference ©• Values for the 
high-mass X-ray binaries (HMXBs) Cen X-3, SMC X-1, and LMC X-4 are from 
Reference {7}, Her X-1 from Reference (JHJ), Vela X-1 from Reference @, and 4U 
1538-52 from Reference ©• Values for the low-mass X-ray binaries (LMXBs) 
Cyg X-2 is from Reference (JTUJ), Cen X-4 from Reference (JTTJ), X1822-371 from 
Reference (JHJ), and XTE J2123-058 from Reference G3J). The two objects plotted 
as crosses instead of dots, the HMXB 4U 1700-37 JHJ) and the LMXB 2S 0921- 
630 ()15|) are blackhole candidates, but may be neutron stars. Error bars are 1<t 
errors. 

limited by the Eddington rate Medd ~ 10 _8 M Q yr _1 (16), the compact stars in 
HMXBs, and consequently also in PSR+NS systems, have a mass very close to 
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Figure 3: The initial mass function of neutron stars as predicted by stellar evo- 
lution theory. The continuous line shows results from Reference (|18[) and the 
dotted line is adapted from Reference (jTTj) . 

their original birth mass, whereas in LMXBs, and consequently in PSR+WD sys- 
tems, the neutron star's mass may have increased significantly. Stellar evolution 
theory Q17|) predicts that stars with a main sequence mass below ~ 18-19Mq 
will produce, when ending their life in a core collapse supernova, proto-neutron 
stars with masses between 1.2 to 1.5M & , whereas more massive ones will produce 
remnants with masses ~ 1.7-1.8M , or a black hole. Fallback during the super- 
nova explosion may alter this (|18|) . In Figure we plot two expected initial mass 
functions of neutron stars. These theoretical theoretical initial mass functions are 
in good agreement with the measurements shown in Figure EJ where PSR+NS 
masses are compatible with all PSR+NS in the range of 1.2-1.5Mq as well as 
HMXBs, with the possible exception of Vela X-l, which may be a representative 
of the second predicted peak at 1.7-1.8 M@. More massive neutron stars are 
expected in LMXBs and their PSR+WD offsprings, which is confirmed by the 
very recent high-mass measurement of PSR J0751+1807 (J3J). 

2.2 Temperature and Radius Measurements 

2.2.1 Elementary modeling of Neutron Star Thermal Emission The 
radiation from a point at the surface of the star is specified completely by the 
specific intensity I(fh,E), defined such that I(rh, E)dCl dE dAdt is the energy 
radiated during a time dt from an area dA as photons of energy in the range 
[E, E + dE], which are emerging in the direction fh within a solid angle dO, ( 16 ). 
If the radiation is polarized, e.g., in presence of a magnetic field, one may sep- 
arate I into two components corresponding to the two polarization modes. By 
integrating over d£l in the outward direction, one obtains the specific, or spec- 
tral, flux F(E) = J J d£l cos 5 1(m, E), where 5 is the angle between m and the 
outward normal n to the stellar surface. Then the flux F = j dEF(E). The 
star's luminosity is then obtained by integration of F over the whole stellar sur- 
face, L = ff dAF = 4ttR 2 F, where F must an average value over the stellar 
surface in case it is not uniform. For blackbody emission and a uniform surface 
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Figure 4: (Left) Effective cross section for interstellar absorption, assuming a 
typical interstellar medium chemical composition. Strong edges from various 
elements are indicated. Figure taken from Reference (|19|). (Right) Observable 
spectrum /oo(-Eoo) from a 1.4 M neutron star with a 10 km radius, i?oo = 13.6 
km, at a distance of 500 pes, assuming blackbody emission at T eoo = 10 6 K. 
Four cases interstellar absorption with Njj = 10 19 , 10 20 , 10 , and 10 22 cm - 2 are 
shown, as well as the case of no absorption. Notice that below 13.6 eV all spectra 
are identical. 

temperature T, we have F = usb?" 4 and hence L = A-kR^u^T^ . For nonblack- 
body or/and a non-uniform emission, it is custumary to express the luminosity 
as L = Air R 2 cfsbT^ , a relation that simply defines the star's effective temperature 
T e . 

Owing to gravitational redshift, a photon emitted at energy E will be detected 
by an observer at infinity with energy Eoo = e®E where (for a non-rotating star 
in the Schwarzschild metric) 

$ _ 1/2 f 2GM 

The luminosity Loo = e 2 *L is redshifted twice because it has both energy and 
time content. Similarly, the fluxes at infinity are F^ = e 4 ®F and F 00 (E 00 ) = 
e 3 *i ? (e~*-E' 00 ). A blackbody at temperature T, once red-shifted, would still be 
seen as a blackbody, but one at temperature = e*T; hence one also defines 
the redshifted effective temperature T eoo = e*T e so that 

Loo = ^Rl F^ = 47T2& ffSB^oo • (2) 

The radius at infinity has to be Rqo = e~®R for consistency, but also has the 
physical interpretation that an observer able to resolve the angular size of the 
star would see it as a disk of radius Roo > R. Observers commonly summarize 
the results of an observation's analysis by giving values of T eoo and Roo- Note 
that a realistic atmosphere spectrum of temperature T e does not redshift into 
the spectrum at temperature T eoo , so that the redshift is in principle measurable 
from spectral fits. 
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Figure 5: (Left) Spectral fluxes of emerging radiation from nonmagnetized at- 
mospheres of various chemical compositions and effective temperatures. (Right) 
Spectral fluxes of emerging radiation from magnetized hydrogen atmospheres of 
various effective temperatures and magnetic field perpendicular to the surface. 
Figure taken from Reference (|2H>. 

Unfortunately, the luminosity , or the flux /oo = L^/iirD 2 = (R^/D) 2 
with D as the star's distance, is not directly observable because of interstellar 
absorption. The probability for a photon of energy E to be absorbed by ionizing 
an atom in the interstellar medium is exp[— Nh<?cR (F)], where Nh is the hydrogen 
column density and a e s(E) is the effective cross section. The quantity that is 
actually observed, the star's spectrum, is hence the specific flux fobs(Foo), given 
by 

Us(Eoo) = e-^^/oc^oo) = e-^^oo) ^V e 3* F(Eooe -*) (3) 

in terms of the specific flux F(E) as emitted by the stellar surface. An example 
of a e ff(E), exhibiting its overall E~ s dependence, is shown in Figure |1J as well 
as the effect of a range of values of Nh on the observable blackbody spectrum 
of a neutron star. Considering that most observed cooling neutron stars have 
temperature of the order of 10 6 K or lower and that a T = 10 6 K blackbody has 
its peak at 0.28 keV, one sees that most photons are absorbed unless the star is 
quite close, i.e. Nh < 10 20 cm~ 2 . 

A spectral fit of / Q b s with a blackbody model involves only three parameters, 
Too, Nh, and Rqo/D, because using the blackbody spectrum Feb for F in Equa- 
tion|3]yields e 3 * Fbb(T; E) = Fbb (Too ; ^oo ) > where e* drops out. For blackbody 
emission and Njj essentially determine the shape of the observed spectrum 
and Roo/D is a scaling factor. Having determined and Nh, if the distance D 
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Figure 6: Expected surface temperature distribution for three different magnetic 
field configurations. In the left panel the field is a pure poloidal dipole with its 
north pole at (9, <fi) = (90°, 90°) and a strength of 10 12 G. In the central panel a 
poloidal quadrupolar component has been added to the poloidal dipole, whereas 
in the right panel a todoidal component, inmersed deep in the crust and with the 
same symmetry axis as the poloidal component, has been added. Figure taken 
from References (|31|) and (|3*2*|). 

is known, one obtains a measurement of Roo. A reliable measurement requires, 
however, that the real spectrum F{E) does not differ too much from that of a 
blackbody. 

2.2.2 Models of The Neutron Star Surface The assumption that the 
locally emitted spectrum F(E) is a blackbody is very dubious. Romani (|2*U)) was 
the first to show that that deviations of F(E) from Fbb(E) can be very large when 
the neutron star surface is an atmosphere. The extent of the deviations depends 
strongly on the chemical composition, which, unfortunately, is unknown. For a 
light-element composition, H or He, the opacity decreases strongly with photon 
energy, ~ E~ 3 , and high-energy photons emerge from deeper, and hence hotter, 
layers: The resulting spectrum has a strong excess compared with a blackbody's 
Wien tail and a strong deficit in the Rayleigh-Jeans regime, as can be seen in the 
examples shown in the left panel of Figure El For an atmosphere composed of 
heavier elements (12211231) , the opacity energy dependence is not as strong and the 
emerging spectrum is closer to a blackbody than in the case of light elements, 
but presents numerous absorption lines, as, e.g., in the Fe atmosphere spectra 
shown in Figure [5J In the presence of a strong magnetic field, as is the case for all 
isolated cooling neutron stars, the specific intensity becomes strongly anisotropic 
(|24j) and the resulting spectra are, usually, intermediate between a blackbody 
and the corresponding nonmagnetized spectra (25,26). Examples of magnetized 
H atmosphere spectra are shown in the right panel of Figure 03 

Given the strong increase in electron binding energies in the presence of a 
strong magnetic field, it is also possible that the neutron star surface may be in a 
condensed state, i.e., a liquid or a solid (f2*7|) . Models of the spectrum emitted by 
such a surface have been presented only recently |28 ), and for some specific field 
strength and chemical composition, it may simulate a blackbody in the observable 
energy range of 0.1-1 keV. 

The presence of a strong magnetic field is not only felt at the surface, but also 
felt much deeper in the star's crust where heat from the interior is transported to 
the surface by electron conduction. Thermal transport transverse to the field is 
strongly suppressed so that heat essentially flows along the magnetic field lines, 
resulting in a nonuniform surface temperature distribution (29,30). When ob- 
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Figure 7: Upper limits on the thermal luminosities of neutron stars in supernova 
remants (names are in parenthesis). Figure adapted from Rerefence (33), with 
Reference (|34|> 

serving the neutron star thermal radiation in the X-ray band, one preferentially 
detects photons from the warm regions of the surface. If extended cold regions 
are present at the surface, they may nevertheless dominate the thermal flux at 
low energies, i.e., in the optical range. Examples of possible surface tempera- 
ture distributions owing to the anisotropic thermal conductivity are illustrated 
in Figure H3 

2.2.3 Temperature and Thermal Luminosity Measurements of Iso- 
lated Cooling Neutron Stars The very brief description presented above 
gives a flavor of what is involved in interpreting observational data about cooling 
neutron stars and what is needed to obtain reliable measurements of their tem- 
peratures, thermal luminosities, and, we hope, radii. We summarize in Figures [7| 
and|Slthe best presently available results. All objects shown as PSR are active 
pulsars and have a very energetic magnetosphere - which is a copious X-ray, and 
in some cases also 7-ray, emitter - producing a power-law spectrum superposed 
on the thermal emission. Separating the surface thermal emission from this mag- 
netospheric emission is not a trivial task. For the Crab pulsar, PSR B0531+21 
and PSR J1124-5916, the magnetospheric emission is so strong that the thermal 
emission is undetectable and only upper limits on the latter are obtainable. In 
the similar case of PSR J0205+6449 in 3C58, the thermal emission is barely de- 
tected as a slight correction to a pure power-law spectral fit. RX J0007. 0+7302 in 
the supernova remnant CTA1 manifests itself as a point source with a power-law 
spectrum embedded in a pulsar wind nebula, and no pulsations have been de- 
tected to date: The thermal emission is also undetected and hence only an upper 
limit on it is possible. The four putative objects in Figure [7| marked "?" come 
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from a deep search for central objects in these four supernova remnants Q35|) : No 
compact object has been found and the quoted upper limit is very restrictive if 
the compact object is a neutron star, but the latter may also be a black hole. 

For the stars in which the thermal spectrum is clearly detected, the dilemma is 
in the choice of the theoretical spectrum to be used in the fits. In the results shown 
in Figure |H] we distinguish between spectral fits performed with blackbody and 
those performed with magnetized hydrogen atmosphere models, depending on 
the resulting Roo estimate. Blackbody fits give higher T eoo and smaller than 
atmosphere model fits: An atmosphere model fit is chosen when a blackbody 
fit gives too small an Roo, and a blackbody fit is chosen when an atmosphere 
model fit gives a much too large Rqq. Moreover, in several of the cases where the 
atmosphere model fit is prefered, the blackbody fit is statistically unacceptable. 
Considering, moreover, that only the hottest part of the star may be detected 
in the X-ray band, the deduced temperature may not be T eoo , as defined in 
Equation [2J The values of we report in Figure |H] take into account this 
ambiguity and are probably more reliable observational values to compare with 
theoretical cooling calculations than T eoo , a point of view we adopt in Section [5J 

The two stars RX J1856. 5-3754 and RX J0720. 4-3125 are two of the Magnificent 
Seven X-ray dim isolated neutron stars, all of which share the peculiarity that 
the blackbody fit of their X-ray thermal spectra results in R^ much smaller than 
10 km (36). However, optical detections show in several cases a clear Rayleigh- 
Jeans spectrum corresponding, when fitted with a blackbody, to a much lower 
temperature than the X-ray-detected spectrum and to a much larger emitting 
area (36,37). They may be candidates for neutron stars having strong internal 
toroidal fields, as illustrated in the right panel of Figure El and reinforce the case 
that, to date, radii inferred from spectra of strongly magnetized neutron stars 
are unreliable. For example, RX J1856. 5-3754 was proposed as a candidate for a 
quark star (|38|). but the optical data do not support an anomalously small radius 

2.2.4 Radius Measurements of Neutron Stars in Quiescent Low- 
Mass X-ray Binaries It is clear from the previous discussion that measure- 
ments of neutron star radii through fits of their thermal spectra will give reliable 
results only if the following three conditions can be met: (a) The chemical compo- 
sition of the atmosphere is known, (b) the magnetic field is small enough, < 10 10 
G so as to not affect the spectrum, and (c) the star's distance can be accurately 
measured. Quiescent LMXBs in globular clusters fulfill these conditions (40). 
Owing to past accretion, H is present and, because of the strong gravity, sedi- 
mentation assures that it will float to the surface, whereas heavier elements will 
sink, resulting in a pure H atmosphere. As in all LMXBs magnetic fields are 
negligibly small and globular clusters, containing ~ 10 5 -10 6 stars, have distances 
that can be measured with an accuracy of 5-10%. A possible drawback is that, 
despite the system being quiescent, it is not yet possible to exclude that accre- 
tion is still occurring at a very low rate, adding constantly to the atmosphere 
heavy-elements that significantly alter the spectrum. Being in themselves very 
interesting systems, many globular clusters have been, and are still observed by 
CHANDRA and XMM-Newton, which found (gU) more than 20 quiescent LMXBs 
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Figure 8: Measured thermal luminosities of isolated cooling neutron stars. The 
names of the associated supernova remnants, when such a remnant is known, are 
given in parentheses. Figure adapted from Reference f3"3"f 

(qLMXBs), and more candidates are constantly being detected. 

The radii of the two qLMXBs in globular clusters u Centauri and M13, which 
show purely thermal spectra, have been recently estimated as = 13.6 ± 0.3 
km (|42l> and Rqo = 12.8 ± 0.44 km (J43J), respectively. However, the fits were 
realized with a library of H atmosphere spectra from Reference 022|) which were 
calculated at fixed surface gravity g s = 2 A3 x 10 14 cm s -2 . Because i?oo and 
e* are parameters of the spectral fit (Equation |3J) and independent variations 
of them are equivalent to independent variations of M and R, a self-consistent 
spectral fit requires a set of model atmosphere spectra F for the whole range of 
fitted g s . Heinke et al. (|44j) recently showed that such a self-consistent analysis 
leads to significant changes in the deduced Roo, and we show in Figure IT51 their 
results for the qLMXB X7 in the globular cluster 46 Tucanae. 

2.3 Transiently Accreting Compact Stars 

The soft X-ray transients (SXRT), also known as X-ray novae (|45|) . are a class 
of LMXBs in which accretion is not continuous. They present repetitive phases 
of high accretion separated by periods of quiescence. During outbursts they 
reach very high luminosities, L Q ~ 10 37 -10 38 erg s , whereas in quiescence the 
luminosity drops by many orders of magnitude, L q < 10 33 erg s . The typical 
duration of an outburst, t Q , is much shorter than the recurrence time between 
outbursts, t r . A common interpretation of this bimodal behavior is the disk 
instability model that was developed to explain dwarf novae 1)46(1 . which are 
similar systems in which the accreting star is a white dwarf, and extended to 
SXRTs Q47() : The companion star is losing mass at a very low rate and feeding a 
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disk that becomes periodically unstable and empties rapidly onto the central star, 
producing a spectacular outburst. Many SXRTs contain black hole candidates, 
but in some cases Type I X-ray bursts occur and we can be certain that the 
accreting object is a compact star. Below we consider only the systems that are 
known to contain a compact star. 

During the quiescent phases between outbursts, many of the SXRTs show a 
thermal X-ray spectrum corresponding to surface temperatures of the order of 
10 6 K: We see the surface of the neutron star that has been heated during the 
accretion phase(s). Given the high internal temperature of the star, the heat 
released in the upper layers by accretion and thermonuclear burning flows back 
to the surface, because of the large temperature gradient in the envelope, and is 
radiated away. Non-equilibrium processes inside the star, induced by the com- 
pression of matter because of accretion, release heat that is stored in the stellar 
interior and that leaks out slowly when accretion stops. Such processes certainly 
occur within the crust: Iron-peak nuclei produced by thermonuclear burning at 
the surface, when pushed to higher densities, undergo a series of reactions - elec- 
tron capture, neutron emission and pycnonuclear fusion Q48|) - until they dissolve 
into the star's core. Overall, an amount of heat, Q nuc ~ 1.5 MeV, is released in 
the crust for each accreted baryon. Brown et al. (|49|) showed that this energy is 
more than sufficient to explain the observed quiescent (thermal) luminosities L q 
in terms of the time-averaged accretion rate, taken over many accretion cycles, 
<M>: 



L q x f Q nuc ^ * 6 x 10 32 / J2=p_ <f> ergs" 1 , ( I) 
m u 1.5 MeV 10 il M yr 1 



where the coefficient / represents the fraction of Q nnc which is actually stored in 
the stellar interior, i.e., not lost by neutrino emission. The luminosity L Q during 
an accretion outburst can be estimated as L ~ (AM/t )(GM/R), where AM is 
the mass accreted during the outburst, and similarly, <M> can be estimated as 
AM/t T . From Eq. 0] one obtains 

i8 . ioX/ ^r%. ioX/ ^A (5 , 

t r V Rc J m u c z 100 to 200 

where the 100 to 200 range corresponds to the possible range of 2GM/Rc 2 from 
~ 0.33 to ~ 0.66, or 

Lq U 

J L (100 to 200) x t y J 

Note that L /L q is independent of the source's distance, which is often poorly 
constrained, whereas t D and t T can in principle be obtained directly by monitoring 
the source for a long enough time. 

We present in Tabled a list of SXRTs for which enough information is available 
to deduce, from Equation H3 an estimate of the storage efficiency /. The six 
systems in the upper part of the Table have been detected in outburst many 
times, so estimates of t Q and t r are reliable, whereas the four cases in the middle 
of the table have exhibited only two or three outbursts and for which estimates 
of / have to be taken with caution. For the five systems in the lower part of 
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Table 1: Storage efficiency of transiently accreting neutron stars 



SoilTVP 


Storage 
efficiency / 


1? pfprpn rp (si 


Aql X-1 


~ 1 


(50) 


CXI (in NGC 6440) 


~ 1 


(ED 


ATT 1 fiDS ^9 




\ v\} ) 




D 1 


( ^9 \ 


EXO 1745-248 


~ 0.01 


JSSJ) 


SAX J1808. 4-3658 


~ 0.01 


(54) 


Cen X-4 


~ 0.1 


(50) 


XTE J1709-267 a 


~ 0.01 


© 


XTE J1751-305 


- 0.1 




XTE J1807-294 


~ 0.1 


(57) 


XTE J0929-314 


t r /780 yrs 


(|56[l57j) 


KS 1731-260 a 


t r / 1,000 yrs 




X 1732-304 a 


t r /500 yrs 




EXO 1747-214 


f r /l,500 yrs 


dnni) 


XTE J2123-058 


i r /67 yrs 





° Quasi-persistent sources: Outbursts last several years instead of only a few weeks to a few 
months. 

the table, only one outburst has been detected, t Q is poorly constrained, and t r 
is at best a guess. Thus we give for / an entry of i r /number-of- years, i.e., a 
best estimate of the right side of Equation |f)J However, for four of these five 
systems, an / of order one would require very long recurrence times t r , which are 
well beyond what the disk instability model can accommodate (|47|). However, 
the validity of the storage efficiency / rests on the assumption that the neutron 
star has reached a quasi-stationary state, which requires fairly regular recurrent 
accretion for at least ~ 10 4 yrs (62 ), a time scale far beyond present coverage of 
the sources ! Notwithstanding these limitations, it is striking that in most cases 
the storage efficiencies seem to be quite low and seem to indicate the occurrence 
of strong neutrino losses. 

2.4 Other Possible Observational Constraints 

The above selection of observational constraints is of course incomplete and re- 
flects our bias in judging relevant and promising probes of the physics of the inner 
core. We very briefly mention here some other possible constraints. 

Quasi-periodic oscillations at millisecond periods observed in several LMXBs 
may constrain the mass and radius of the neutron star, but interpretation of 
the origin and nature of these phenomena is still controversial and inconclusive 
(|63|). Detection of atomic absorption lines in the thermal spectrum of a neutron 
star can directly yield the gravitational redshift due to the star when the line 
is identifiable. The identification can be secure only for stars with a magnetic 
field sufficiently weak so that it does not affect the atomic structure. Candidates 
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for these measurements are the neutron stars in LMXBs, and there is, to date, 
only one case of such detection: Fe XXVI and XXV for n=2-3 and O VIII for 
n=l-2 transitions in the burst spectrum of EXO 0748-676 (64), which results in a 
redshift of z = 0.35, i.e., e~* = 1.35 1 . Also, in Type I X-ray bursts from LMXBs, 
fits of model spectra to the observed burst spectra can potentially measure both 
M and R (65 ), but model spectra do not seem to be accurate enough to produce 
reliable and reproducible results. The recently discovered superbursts (|6l))) are 
interpreted as the step following the normal short Type I X-ray bursts, which are 
due to explosive burning of 4 He into 12 C at the surface of an accreting neutron star 
in a LMXB, which consists of explosive thermonuclear burning of the accumulated 
layer of C into iron-peak nuclei. However, reproducing the observed properties of 
these superbursts imposes very tight constraints on the temperature in the crust 
of these stars and on the neutrino emission from the core Q67|) . 

In isolated neutron stars, two more constraints emerge. From the observations 
of glitches in several young radio pulsars, Link et al (6SJ) have deduced, under the 
commonly accepted scenario that these involve differential rotation of the neutron 
superfluid in the inner crust, that the relative moment of inertia of the crust is 
at least 1.4% of the total stellar moment of inertia. The second constraint comes 
from the recently discovered pulsar PSR J1748-2446ad (|69|) . the fastest known 
pulsar to date with a period of 1.39 ms. For a given assumed mass, it imposes the 
stringent limit that its equatorial radius must be smaller than the mass-shedding 
radius, at which point the centrifugal force becomes larger than gravity. Both 
results translate into M versus R constraints (6 8Tl70f) . which we use in Section [23 
and Figure IT51 

3 DENSE NUCLEON MATTER 

Over the past decade there have been numerous attempts to compute the bulk 
properties of nuclear and neutron-rich matter. These include microscopic many- 
body calculations using realistic nucleon-nucleon potentials and phenomenologi- 
cal relativistic and non-relativistic mean-field theories. The former approach em- 
ploys a potential that provides an accurate description of the measured nucleon- 
nucleon scattering data and uses variational or quantum Monte Carlo techniques 
to obtain the EoS (for a review of these methods, see Reference (|71|). and for 
recent progress in using variational techniques to obtain the EoS of relevance 
to neutron stars, see Reference (|72|)). Variational techniques rely on the pa- 
rameterization the nuclear wave-function. Correlations induced by interactions 
between pairs of nucleons are incorporated through suitably parameterized pair- 
correlation operators, which act on the symmetrized product of the Fermi gas 
wave-function. In contrast, quantum Monte Carlo treatments such as the Green's 
function Monte Carlo are not limited by the form of the variational wave- function 
and can potentially include all possible correlations in the many-body system. 
This approach is, however, computationally intensive owing to the spin and 
isospin dependence of the nuclear interaction and grows roug hly as 2 A A\/(N\Z\) 
where A, N, and Z are the baryon number, neutron number, and proton number, 

*Fe XXV and Fe XXVI are Fe with only two and one electrons bound to it, respectively 
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respectively. Although this method suffers, in principle, from the fermion sign 
problem, Carlson et al. (|73j) have developed algorithms to minimize its influence 
and have calculated the ground-state energy of 14 neutrons in a periodic box. 

The other notable method that employs realistic nucleon-nucleon potentials 
is the Brueckner-Hartree-Fock (BHF) approach. Here the bare nucleon-nucleon 
interaction is used to determine the interaction energy for pairs of nucleons. The 
pairs are treated as independent particles, and correlations that arise because of 
Fermi statistics are incorporated. This method, developed by Brueckner, Bethe, 
and Goldstone, is nonperturbative in the coupling but utilizes a perturbative ex- 
pansion in the number of independent hole lines to make it tractable (|74|) . Recent 
studies have shown that the convergence of this expansion can be improved if the 
single-particle dispersion relation is obtained using the Hartree-Fock approxima- 
tion. In the lowest-order Brueckner-Hartree-Fock (LOBHF) method, only the 
lowest-order terms in the hole-line expansion are retained. The LOBHF results 
are in fair agreement with those obtained using variational methods (75). 

The realistic nucleon-nucleon potentials are well constrained by a collection 
of scattering data compiled by the Nijmegen group (J7BJ). However, these poten- 
tials are not unique. The data probes the scattering for energies below the pion 
threshold of 350 MeV. Thus, various models for the nucleon-nucleon interaction 
with different short-distance behavior can be constructed to reproduce the elas- 
tic scattering phase shifts up to the measured energy. These different potentials 
are equivalent from the point of view of the low-energy data. Alternatively, low- 
energy observables cannot be sensitive to the details of short-distance physics. 
Nonetheless, these potentials can differ in their predictions for the many-body 
system because additional effective three-, four-, five-, etc. body forces can be 
relevant. The strength of these many-body forces is not unique and depends on 
the assumed form of the two-body force at short distances. Fortunately, for the 
realistic two-body potentials employed to date, the addition of three-body forces 
is sufficient to describe light nuclei Q77j) . Several groups developed an alternate 
description of nucleon-nucleon interactions based on effective field theory (for a 
review, see Reference 1)781)). This approach allows for a systematic treatment of 
the short-distance physics by organizing the calculation in powers of momenta 
and identifies a small expansion parameter at small momenta. It obviates the 
need to make specific assumptions about the short-distance potential. Bogner et 
al. ( 79 ) have shown that a momentum-space potential with a rather small cutoff 
of ~ 2 fm _1 can describe low-energy data. Furthermore, they have shown that all 
realistic models of the nucleon-nucleon potential evolve (in the renormalization 
group sense) to this form. Although it is not clear if these low- momentum interac- 
tions will make many-body calculations more tractable (or less non-perturbative), 
preliminary work is promising ( 80 ) . 

In the traditional nuclear mean-field models, the relation to nucleon-nucleon 
scattering is abandoned in favor of a phenomenological interaction whose parame- 
ters are determined by fitting the model predictions (in the mean-field approxima- 
tion) to empirical properties of bulk nuclear matter at nuclear saturation density. 
These mean-field models, such as the non-relativisitic Skyrme model and the 
relativistic Walecka model and its variants, can be viewed as approximate imple- 
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mentations of the Kohn-Sham density functional theory, which has been widely 
studied and used in quantum chemistry and condensed matter physics (81). 

Empirical properties of bulk nuclear matter are extracted from the analysis 
of nuclear masses, radii, and excitations of large nuclei. These analyses reliably 
disentangle the nuclear bulk properties from surface and electromagnetic con- 
tributions. Without delving into the details and caveats of these analyses, we 
simply state the empirical properties for which there is broad consensus: (a) 
Nuclear saturation density (the density at which symmetric nuclear matter is 
bound at pressure P = 0) is po = 0.15-0.16 nucleons fm -3 ; (b) the binding 
energy per nucleon at saturation is 16 MeV; (c) the nuclear compression mod- 
ulus defined through the relation Kq = 9dP/dp\ p=Po is determined to be in the 
range of 200-300 MeV; (d) the nucleon effective mass at saturation density is 
0.7 — 0.8 ttt-tv; and (e) the nuclear symmetry energy defined through the relation 

5 = 1/2 d 2 (e/ p)/d(l — 2x p ), where x p is the proton fraction, is in the range of 
30-35 MeV. This empirical knowledge provides valuable constraints and guidance 
for models of dense nuclear and neutron-rich matter. 

3.1 A Simple Model for Nuclear Matter 

To provide a description of nuclei and nuclear matter, Walecka (|82|) proposed a 
field-theoretical model in which the nucleons interact via exchange of scalar and 
vector mesons. The model has been refined and used extensively to study nuclear 
properties in the mean- field approximation (for a review, see Reference (83)). 
With the intent of providing a pedagogic overview of the various forces at play in 
dense nuclear matter, we describe briefly one such model introduced by Boguta 

6 Bodmer (|51)) (a more elaborate discussion of this model and its application to 
the structure of neutron stars can be found in Reference ( 85)). 

The model proposes that in dense matter, nucleons interact with effective short- 
range forces. The Lagrangian is given by 

Cn = ~*n (iVfy ~m* N - g^Nl^ ~ g P Nl^r N ■ R^ V N 

+ \d^o - \mlo 2 - U(a) - 

+ \mlV^ - X -R pv ■ RT + \m 2 p R p ■ fr, (7) 

where m* N = — g a NO is the nucleon effective mass, which is reduced in 
comparison to the free nucleon mass owing to the scalar field a, taken to 
have m a = 600 MeV. The vector fields corresponding to the isoscalar omega 
and isovector rho mesons are given by V pu = d p V u — d u Vfj,, and R pv = d p R v — 
d v R p respectively. The exchange of these mesons mimics the short-range forces 
between nucleons. In addition to the coupling between nucleons and mesons, a 
self-interaction between scalar mesons given by 

b c 
U(a) = -m N (g cjN af + -(g aN a) 4 , (8) 

where b and c are dimensionless couplings, is introduced to obtain good agreement 
with the empirical value of the nuclear compressibility (|84[). is the nucleon 




Figure 9: The nuclear equation of state. (Left) The density dependence of the 
baryon chemical potential, the electron chemical potential, and the energy per 
baryon. (Right) The relation between energy density and pressure. 



field operator and tn is the nucleon isospin operator. The five coupling constants, 
9a N, 9uN, 9pN, b, and c, are chosen, as in Reference (|85]> . to reproduce five 
empirical properties of nuclear matter at the saturation density listed above. 

The model is solved in the mean-field approximation. Here, only the time 
component of the meson fields have nonzero expectation values. The symbols a, ui 
and p denote sigma-, omega-, and rho-meson expectation values that minimize 
the free energy given by (|55]l 

^nuclear On,/"e) = ij dk k 2 [e n (k) - /i n ] + J dk k 2 [e p (k) - fip] 

+ l(mla 2 -m 2 ^ 2 -m 2 pP 2 )+U(a)-J^, (9) 



where 



= \jk 2 + m* N 2 + g wN uj - ^g p NP , (10) 



e p (k) = \Jk 2 + m* N 2 + g uN u + ^g p NP , (H) 

are the neutron and proton single-particle energies. The single-particle energy 
at the Fermi surface defines the respective chemical potentials. In stellar mat- 
ter these chemical potentials are related by the condition of weak interaction 
equilibrium (n + v e <-> e~ + p) given by 

Pn= P-p + P-e (12) 

(note that we have set p Ue = since neutrinos are not typically trapped). Con- 
sequently, only pL n and p, e are independent and correspond to the two conserved 
charges, namely baryon number and electric charge. Furthermore, we require 
that bulk matter be electrically neutral. To enforce local charge neutrality, we 
require the charge density pq = dQ QUC \ eajI /dp e = 0, which uniquely determines 
p e . Therefore, in effect only p n is independent and the free energy or pressure as 
a function of p n completely specifies the EoS of dense matter. 
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Figure El shows the thermodynamic properties of charge- neutral stellar matter. 
The left panel shows the baryon chemical potential, electron chemical potential 
and the energy per particle as a function of the baryon density. The right panel 
shows the energy density as a function of the pressure and is usually referred 
to as the equation of state (EoS). This relation that is required to solve for the 
structure of the neutron star and determines the neutron star mass and radius. 
An EoS that has on average a larger (lower) pressure for a given range of energy 
density is termed stiff (soft). 

Different nuclear EoSs constructed to satisfy the empirical constraints at nu- 
clear density and hence similar in this regime can differ significantly at lower and 
higher densities. In addition, the difference between symmetric nuclear matter 
EoS and the neutron-rich stellar matter EoS could be significant. This difference 
arises mainly because of the difference in the density dependence of the nuclear 
symmetry energy. The magnitude and the density dependence of the proton frac- 
tion in particular is sensitive to it. In the mean-field model considered here, the 
nuclear symmetry energy arises owing to the isovector force from the exchange of 
p mesons, and its density dependence is linear. In contrast, in more sophisticated 
treatments that employ realistic nucleon-nucleon interactions and correlations be- 
yond mean-field theory, the symmetry energy has non-trivial density dependence 
and its magnitude is typically smaller. 

These differences lead to two generic trends that distinguish variational treat- 
ments such as those reported in the work of Akmal, Pandharipande, and Raven- 
hall (APR) (O from the mean-field EoSs. First, the APR EoS for beta-stable 
matter is considerably softer than the mean-field EoS at low density and stiffer 
at high density. Second, the typical proton fraction at high density in the APR 
EoS is smaller than those predicted by the mean- field EoS. Both of these trends 
lead to important consequences for neutron stars. The differences in the low- and 
high-density behavior of the EoS, as we discuss in more detail in Section lead 
to more compact (smaller radii) and more massive (larger maximum mass) neu- 
tron stars in the APR EoS. Because these differences can in part be attributed to 
differences in the density dependence of the symmetry energy, ascertaining this 
directly from experiment would be desirable. Recent work has emphasized that 
this can be extracted from precise measurements of the neutron skin of heavy 
nuclei (|5B|) and heavy- ion collisions (|57j). 

The lower proton fraction seen in APR suppresses neutrino cooling of the 
core. This occurs because the single-nucleon neutrino-producing reactions such 
as n — > p + e~ + i7 e and e~ + p — > n + v e , which are termed direct URCA 
reactions, are kinematically forbidden when the proton fraction is small. These 
reactions, which occur on the Fermi surface, can conserve momentum only if 
\kFn\ < |&Fp + kFe\ (termed the triangle inequality). Because electric neutrality 
requires that |A:_F e | ~ \kF P \, it follows that the triangle inequality is difficult to 
satisfy when the proton fraction is small. A detailed calculation shows that the 
direct URCA reaction occurs when the proton fraction is x p = n p j (n n +n p ) > 1/9 
(|88|) . For smaller values of x p , the dominant neutrino-producing reactions involve 
two nucleons in the initial and final state to ensure momentum conservation. They 
include the charged-current processes nn — > npe~v e and np — ► ppe~D e (and their 
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inverse), which are termed modified URCA reactions, and the neutral-current 
processes such as nn — > nnvu, which are termed neutrino bremsstrahlung. These 
two nucleon reactions rates are suppressed by a factor of (T/Tp) 2 in degenerate 
matter relative to the direct URCA reaction. They are roughly five orders of 
magnitude slower for the typical temperatures of interest in neutron star cooling 
(see Section 5). 

3.2 Nucleon Superfluidity and Superconductivity 

The low-energy excitation spectrum and the response properties of the dense 
nuclear matter play key roles in the dynamics and thermal evolution of compact 
stars. As we discuss in Section |21 young neutron stars cool primarily owing 
to neutrino emission from the core. Thus, phenomena such as superconductivity 
and superfluidity, which drastically alter the low-lying spectrum, become relevant 
even if their influence on the EoS is not significant. Much like electron pairs, via 
the Bardeen-Cooper-Schreiffer (BCS) (89) mechanism in low-temperature solids, 
nucleons can pair because of their intrinsic attractive interactions in nuclei and 
in dense matter. The typical nuclear pairing energy, or gap, is of the order of 
an MeV. Because the ambient temperature inside neutron stars is <C MeV at all 
times subsequent to a few tens of seconds after their birth, nucleon pairing can 
dominate the low-energy energy response properties of the neutron star interior. 

In the neutron star crust, as nuclei become increasingly neutron rich, neutrons 
drip out of nuclei at a density of pdn P — 4.3 x 10 11 g cm~ 3 . At these low densities 
the neutron-neutron interaction is attractive in the 1 So channel, leading to a large 
scattering length of a ~ 20 fm. This attractive interaction destabilizes the Fermi 
surface and results in pairing and superfluidity. The pairing gap calculated using 
BCS theory is given by 

A= t exp Uo)W r ) )' <i3) 

where M is the neutron mass, N(0) = M hp /it the density of states at the Fermi 
surface and V{kp) the interaction potential evaluated at the scale of the Fermi 
momentum. Clearly, the BCS formula provides only an estimate and cannot 
reliably calculate the magnitude of the gap in the strongly interacting system. 
Nonetheless, Eauation ll3l illustrates several key trends seen in more sophisticated 
treatments (for a recent review of superfluidity in dense nucleon matter, see Ref- 
erence (|90j0 . First, at low density where the range of the interaction is unim- 
portant and only the scattering length is relevant, we can write V(kp) ~ ixajM. 
Inserting this into Equation El we see that at low density the gap increases with 
increasing density. At higher density the range of the interaction becomes rele- 
vant, the nucleon-nucleon interaction changes character and eventually becomes 
repulsive. Consequently, the gap will eventually decrease in magnitude and van- 
ish. This trend results in a bell-shaped curve for the gap as a function of density. 
Detailed calculations discussed in Reference (90) indicate that the neutron 1 So 
gap reaches its maximum value ~ 1 MeV when the neutron Fermi momentum 
k-Fn — 0.8 MeV, and vanishes when /cp n > 1.5 fm. 
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Deeper, inside the core, So proton superconductivity and neutron pairing in 
the 3 P2 — 3 -F2 channel become possible. The typical behavior of the proton super- 
fluid gap and its dependence on the proton Fermi momentum is similar to that of 
1 Sq neutron pairing because the interaction and the momentum probed are nearly 
identical. However, because the protons coexist and interact strongly with the 
dense neutron liquid, this changes their dispersion relation and consequently the 
density of states at the Fermi surface. In the simple mean-field models, this is 
encoded in the proton effective mass which is significantly reduced. Typically, 
this leads to a factor of 2-3 reduction in the maximum value of 1 So proton gap 
relative to the 1 So neutron gap. The situation with the 3 i 3 2 — 3 -^2 neutron pairing 
in the core is much less understood. Even simple BCS estimates that ignore any 
form of medium-polarization effects, such as induced interactions and particle- 
hole renormalization of the bare potential, lead to vastly different predictions for 
this gap at high density because modern models of the nucleon-nucleon potential 
are not constrained by data at these larger relative momenta. This is further ex- 
acerbated when medium-polarization effects are taken into account. Earlier work 
indicated that medium-polarization could strongly enhance the gap (|91|) , whereas 
a recent calculation indicates that a medium-induced spin-orbit interaction could 
lead to a large suppression (92). 

BCS pairing between nucleons leads to an energy gap A in the excitation 
spectrum. Consequently, neutrino processes are exponentially suppressed when 
T -C A owing to the paucity of thermally excited quasi-particles. This would 
lead us to naively conclude that nucleon superfluidity leads to an exponential 
suppression of the neutrino cooling rates in the core, where neutron and proton 
pairing is likely. However, in the vicinity of the critical temperature, the system 
is characterized by strong fluctuations as Cooper pairs form and break. These 
fluctuations give rise to a very efficient neutrino-emission process termed the pair 
breaking and formation, or PBF, process (I93U94[) . In the BCS approximation for 
the response function, the PBF process emissivity turns on just below T c , grows 
as T decreases, and eventually becomes exponentially suppressed when T <C T c , 
For typical values of the P2 gap, one finds that at T ~ 10 9 K, the PBF rate can 
be up to one order of magnitude higher than the modified URCA rate (33). 

Many-body correlations beyond those induced by the pairing play a role in 
determining the response properties, especially the weak interaction rates. The 
effects of these correlations on the neutrino emissivity are poorly known. For 
example, until recently the modified URCA reaction and neutrino bremsstrahlung 
reactions were calculated using a simplified one-pion exchange interaction in the 
Born approximation (J95|). This calculation has been revisited using soft radiation 
theorems and a more realistic and nonperturbative treatment of the nucleon- 
nucleon interactions (|96 D . These results indicate a suppression by a factor of 2-3 
in the neutrino rates relative to the calculation in Reference (95). The effects 
due to medium polarization have also been investigated, but the findings have 
been inconclusive. Researchers suspect that medium-polarization effects through 
particle-hole correlations lead to a softening of the pion dispersion relation and 
are relevant for neutrino emissivity (98). In a model calculation, this softening has 
enhanced the modified URCA rate at nuclear density by an order of magnitude 
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(|97|). However, there is still no consensus about the role of these correlations. For 
example, in Ref. (|99|) the authors find that particle-hole screening of the nuclear 
interaction leads to a suppression in the neutrino bremsstrahlung rate ( 99 ) . 

From the preceding discussion it should be clear that there are many open issues 
relating to nucleon superfluidity in the core. The importance of the PBF pro- 
cess, which operates at relatively small values of the gap and dominates neutrino 
cooling, emphasizes the need to know the precise behavior of the 3 i"2- n eutron 
and 1 5o-proton pairing deep inside the core. It is also important to revaluate the 
role of correlations beyond BCS mean-field theory, both on the low-temperature 
specific heat and the neutrino emissivities. 

4 NOVEL PHASES 

With increasing density the chemical potential for baryon number and negative 
electric charge increase rapidly owing to the repulsive nature of strong interactions 
between nucleons at short distances. This furnishes energy for the production of 
strange baryons and the condensation of mesons. At higher densities, our knowl- 
edge of quantum chromodynamics (QCD) and its asymptotic behavior leads to 
the expectation that quarks inside nucleons will delocalize and form a uniform 
Fermi sea of quarks through a deconfinement transition. These expectations are 
borne from model descriptions. Unfortunately, these model predictions are diffi- 
cult to quantify and the precise density at which these transitions may occur is 
poorly understood. Given the current state of the art, it appears that theoretical 
advances alone are unlikely to provide conclusive evidence for a phase transi- 
tion at supranuclear density. Ultimately, only through their possible presence in 
compact stars can hope to learn of their existence. This would be possible if we 
could (a) elucidate in a model-independent manner generic properties of these 
novel phases and (b) identify and calculate some properties that would affect 
observable aspects of the compact star structure and evolution so that they are 
distinguishable from those predicted for dense neutron-rich matter. Below we 
consider a few of these scenarios and discuss key bulk and response properties 
that distinguish these phases. 

4.1 Hyperons 

Fig. I1UI shows the variation of the chemical potential associated with neutral and 
charged baryons. The thicker curves are predictions of the mean-field model de- 
scribed in the previous section and the thin lines correspond to the case wherein 
strong interactions are ignored. The electrons in dense nuclear matter provide 
a source for the charge chemical potential. Thus the chemical potential for a 
charged baryon B T is n(B T ) = fig ± fi e , where [is = Mn is the baryonic chemical 
potential and /x n is the neutron chemical potential. Consequently the source for 
baryons with different charge differs and negatively (positively) charged baryons 
are favored the most (least). The horizontal dashed lines indicate the vacuum 
masses of the A, S 1 * 1 hyperons. The point at which the corresponding chemi- 
cal potentials cross the vacuum masses are also indicated both for the interact- 
ing and non-interacting nuclear cases. These points correspond to second-order 
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phase transitions at which specific hyperon species begin to be populated. If the 
strong interactions between hyperons and nucleons were attractive (repulsive) 
the second-order transition would occur at lower (higher) density. It is also clear 
from the figure that strong interactions between nucleons play an important role 
as it shows that for non-interacting nucleons the transitions occur at significantly 
larger density. Hyperon-hyperon interactions can also become relevant especially 
if they are strongly attractive as this could lead to a first-order transition at lower 
density (fTUU|) . 




Figure 10: Baryon chemical potentials in dense stellar matter. 



To infer if hyperons exist inside neutron stars it is crucial to either measure 
or compute directly from QCD the hyperon-nucleon potential. Experimental 
data on hyperon-nucleon interactions are very scarce. However, recent develop- 
ments in studying baryon-baryon interactions using lattice QCD are promising 
and could potentially provide useful information in constructing a low-energy ef- 
fective theory for hyperon-nucleon and hyperon-hyperon interactions (|101|) . Ref. 
(102) analyzed various experimental constraints on the hyperon-hyperon interac- 
tion and Ref. (|103|) studied constraints arising from the binding energy of the A 
in hypernuclei. This information was first employed in the mean-field models by 
Glendenning (|104j) . Subsequently a more detailed analysis using models for the 
hyperon-nucleon potentials which are consistent with A binding energy and the 
hyperon scattering lengths have been employed in relativistic BHF calculations 
(105). These studies indicate that the lightest hyperons could appear at relatively 
low density (2-3 times nuclear density). 

The presence of hyperons typically softens the equation of state and enhances 
the response of dense matter. This occurs because relative to the nucleons the 
hyperons contribute more to the energy density than to the pressure because 
they have larger masses and smaller Fermi momenta. They furnish new degrees 



24 



Page & Reddy 



of freedom (that are less degenerate than the nucleons) and can readily partic- 
ipate in excitation processes and thereby enhance the response. In particular 
their presence enhances the neutrino cooling in the core because: (1) they can 
participate in rapid direct URCA reactions such as A — > pe~v e (not kinematically 
suppressed) ; and (2) they reduce the neutron to proton ratio and facilitate the 
momentum conservation needed to initiate the direct URCA involving nucleons 

4.2 Kaon Condensation 

A large number of electrons are required to ensure charge neutrality in dense 
nuclear matter. The typical electron chemical potential /J, e ~ 100 MeV at nuclear 
density. As density increases fi e increases to keep pace with the increasing proton 
number density. The magnitude of this increase is sensitive to the density depen- 
dence of the isovector interaction contribution to the nuclear symmetry energy. 
In Figure ^2 the density dependence of the electron chemical potential in the 
mean-field model is shown as the thicker curve. For reference, the electron chem- 
ical potential for the case of noninteracting nucleons is also shown (thin curve). 
If the energy of a zero momentum negatively charged boson in the medium is 
less than the electron chemical potential it will condense. The amplitude of the 
condensation will in general be regulated by the repulsive interactions between 
bosons. 

In the hadronic phase the likely candidates for condensation are ir~ and K~ . 
In vacuum, pions are significantly lighter than kaons but this situation may be 
reversed in the dense medium owing to strong interactions between mesons and 
nucleons. The physical basis for this expectation is that the effective theory of 
meson-nucleon interactions must incorporate the repulsive s-wave interactions 
arising owing to the Pauli principle between nucleons and mesons composed of 
only up and down quarks, whereas interactions between mesons with u or d 
quarks and strange quarks must be attractive. Experiments with kaonic atoms 
lend strong support to the aforementioned theoretical expectation (jl07[) . In Fig- 
ure ^2 the vacuum pion and kaon masses are shown as the dark and light dashed 
lines respectively. If the masses do not change in the medium, the figure indi- 
cates that 7T~ condensation occurs in the vicinity of nuclear density whereas K~ 
condensation does not occur for the densities considered. When interactions with 
medium are included, a uniform charged pion condensate is disfavored owing to 
a weak and repulsive s-wave interaction. Instead, a spatially varying condensate 
can be favored because of attractive p-wave interactions (for a review see Ref. 
(108)). The kaon-nucleon interaction, however, is strongly attractive. For this 
simple reason we discuss kaon condensation below. 

Kaplan and Nelson p09| proposed the idea that kaons could condense in dense 
nuclear matter. Using a simplified SU(3)x SU(3) chiral Lagrangian they showed 
that K~ could condense at a density approximately three times nuclear density. 
Subsequently, several authors have studied in detail the nature and the role of 
kaon condensation in neutron star matter (for a recent review see Ref. piOj) ). 
Here, we employ a simple, schematic, potential model for kaon-nucleon interac- 
tions considered in Ref. 1)111(1 to illustrate the salient features. The scattering 
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Figure 11: Electron chemical potential, the pion and kaon vacuum and in- medium 
effective masses in dense stellar matter 



length aj(- n characterizes the low-energy kaon-nucleon interaction and experi- 
ment indicates that a^- n ~ —0.4 fm. Following Reference Qlllj) . the interaction 
potential can be modeled as a square well. The parameters of the potential, i.e., 
the depth Vq = —122 MeV and range R = 0.7 fm, are determined by fitting 
the scattering predictions to low-energy nucleon-kaon data (see Ref. ([111[> for 
further details). Given this potential the change in the effective mass of K~ in a 
low-density medium of neutrons can be computed using the Lenz approximation. 
In this approximation, the attractive potential energy experienced by a kaon at 
rest can be related related to the scattering length and is given by 

Vhenz = a K - n p n , (14) 

A* 

where p, is the reduced mass of the neutron-kaon system and p n is the neutron 
density. The effective mass of the kaon computed using Equation El is shown 
in Figure 1111 In this case the kaon effective mass equals the electron chemical 
when p ~ 3 — 4 times nuclear density. This corresponds to the critical density for 
kaon condensation. At higher densities, the Hartree or mean-field approximation 
should be valid. In this case the attractive potential energy of the kaon cannot be 
related directly to on-shell low-energy kaon-nucleon scattering data. The Hartree 
potential is given by 

Air o 

^Hartree = — V R p n , (15) 

where Vo and R are the depth and range of the nK~ potential, respectively. From 
Figure ^2 we see that in the Hartree approximation kaons would condense when 
p ~ 5 - 6 po ■ 
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The presence of kaons inside neutron stars influences their structure and evolu- 
tion in a manner similar to that of hyperons. The softening of the EoS owing to 
a kaon condensate is typically larger than the softening by hypernos because the 
zero-momentum kaon condensate contributes to the energy density but not to the 
pressure. By furnishing negative charge it favors a more isospin-symmetric nu- 
clear phase containing nearly equal numbers of neutron and protons which results 
in additional softening. As we discuss in ^SJ a strong softening in the high-density 
EoS will result in a lower maximum mass and a significantly smaller canonical 
radius. The presence of either pion or kaon condensates result in enhanced neu- 
trino cooling In the case of a kaon condensation enhanced cooling occurs 
owing to both reactions involving kaon decays in the presence of a bystander nu- 
cleon and the nucleon direct URCA reaction made possible by the larger proton 
fractions present in the kaon-condensed phase (|113f) . 

4.3 Normal and Superconducting Quark Matter 

The occurrence of novel hadronic phases depends on the nature of hadronic inter- 
actions and their many-body descriptions. In contrast, the asymptotic behavior 
of QCD, which requires that interactions between quarks become weak with in- 
creasing momenta, strongly supports that at sufficiently high densities nucleonic 
degrees of freedom must dissolve to form a gas of weakly interacting quarks. The 
precise location of this phase transition will depend on model descriptions of both 
the nuclear and quark phases. This is because all model studies indicate that the 
phase transition occurs at rather low densities where perturbative methods do 
not apply. The bag model provides a simple description of the quark matter 
and the confinement (|114j) . The model was designed to provide a description of 
the hadron mass spectrum. The basic tenants of the model are a non-trivial vac- 
uum and nearly free-quark propagation in spaces (Bags) wherein the perturbative 
vacuum has been restored. This restoration costs energy because it requires the 
expulsion of the vacuum condensates. The restoration energy per unit volume is 
termed the bag constant and is denoted as B. The model also provides a very 
simple and intuitive description of bulk quark matter. The pressure in the bulk 
quark phase containing up (u), down (d) and strange (s) quarks is due to the ki- 
netic energy density of quarks and a negative Bag pressure. At zero temperature 
this is given by 

PBa,g(fJ>u, Vd, Vs) = ~ 22 / 7^)3 (V^ 2 +m i ~ B > ( 16 ) 

i=u,d,s ^ 

where 7 = 2 sp ; ns x 3 co ior is the degeneracy factor for each flavor, and /ij and kfi 
are the chemical potential and Fermi momentum of each flavor, respectively. In 
the limit of massless quarks and a common chemical potential for all the quarks, 
the pressure has the simple form 

PbM = ^~B. (17) 

The pressure of bulk quark matter computed using Equation d for B 1 ^ = 150 
MeV and B 1 ^ = 200 MeV are shown in Figure 1121 The leading-order effects 




of perturbative interactions between quarks can also be incorporated in the Bag 
model at high density. This has the effect of renormalizing the kinetic term. At 
densities of relevance to neutron stars, perturbative expansion in a s = g 2 /Air, 
where g is the QCD coupling constant, is not valid. Note that the order a 2 
calculation of the free energy predicts a behavior similar to that of the bag model 
(|115U116|) . Recently, Fraga, et al. recomputed the equation of state of massless 
quark matter to 0(a 2 ). They found the perturbative result is well reproduced by 
the following, bag model-inspired, form for pressure over a wide range of densities 
relevant to neutron stars (|116[) 

3 

^perturb(0(c^), Ll) = — a eff /J 4 - B eS (18) 
47T 

where a e g = 0.628 and B 1 ^ = 199 MeV with a specific choice for the renor- 
malization subtraction point, A = 2/x. The pressure obtained using Equation 1181 
is shown in Figure 1121 For reference, the nuclear APR EoS is shown and the 
crossing points denote the locations of the first-order transitions between nuclear 
matter and quark matter. In the extreme case of the bag model with B 1 ^ = 150 
MeV we see that quark matter is preferred over nuclear matter everywhere. In 
this scenario, quark matter is the true ground state of matter even at zero pres- 
sure and is the basis for the strange quark matter (SQM) and strange quark star 
hypotheses which we discuss below. 
4.3.1 Color Superconductivity. Since the early work of Bardeen, et 
al., it has been well known that degenerate Fermi systems are unstable in the 
presence of arbitrarily weak attractive interactions at the Fermi surface ( 89 ) . The 
instability is resolved by the formation of a Bose condensate of Cooper pairs. For 
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the case of charged fermions, such as electrons, this leads to superconductivity. 
By analogy, the presence of attractive interactions between quarks will lead to 
pairing and color superconductivity. Barrios and Frautschi (|117|) realized this 
several decades ago. The recent realization that the typical superconducting 
gaps in quark matter are larger than those predicted in these earlier works has 
generated renewed interest. Model estimates of the gap at densities of relevance 
to neutron stars suggest that A ~ 100 MeV when \i = 400 MeV (jHHj) . Here, we 
provide a brief introduction to the subject and emphasize aspects that will impact 
neutron star phenomenology (see Reference H119|) for a more comprehensive recent 
review) . 

One-gluon exchange between quarks is attractive in the color antisymmetric 
(triplet) channel and repulsive in the color symmetric (sextet) channel. The 
attraction in the triplet channel can result in s-wave pairing between quarks 
in spin-zero and spin-one channels. Explicit calculations show that the pairing 
energy is significantly large for the spin zero case. This type of pairing occurs only 
between dissimilar flavors of quarks to ensure that the pair (which is a boson) 
wave-function is symmetric. For three massless flavors the condensation pattern 
that minimizes the free energy is known as the color-flavor- locked scheme (120). 
The non-zero condensates in this phase are given by 

(Cl<l^) = -(^<R^) = Ae^e a , M (19) 

where a and (3 are color indices, a and b are flavor indices, and i and j are spinor 
indices. Because e a,/3 ' A e a ^,A = S^d^ — 5£5a it follows that condensation locks 
color and flavor indices. 

In the color-flavor-locked phase all nine quarks participate in pairing. Conse- 
quently there is an energy gap ~ 2A in the fermion excitation spectrum. The 
gluons acquire a mass via the Higgs mechanism by coupling to the colored conden- 
sates. The lowest excitations in this phase correspond to the Goldstone bosons. 
The number and spectrum of Goldstone bosons can be understood by noting 
that the condensate breaks baryon number and chiral symmetries. Chiral sym- 
metry is broken by a novel mechanism |120) because independent Left-handed 
(L) and right-handed (R) condensates are coupled by the vector interaction. The 
octet of flavor Goldstone bosons that result from the breaking od chiral symme- 
try acquire a mass because quark masses explicitly break chiral symmetry. The 
quantum numbers of these pseudo-Goldstone modes map onto the meson octet 
in vacuum. For this reason they are commonly referred to as the pions and kaons 
(or mesons) of the high density phase. Their masses have been computed in ear- 
lier work by matching the effective theory to perturbative QCD at high-density 
(|121[) . Unlike in vacuum, the square of meson masses was proportional to the 
product of quark masses . This resulted in an inverted hierarchy, in which the 
pions were heavier than the kaons. Explicitly, at high density the masses are 
given by 

mt± = a(m v + mj)m, 

4 ± = n L:+4«: 

m 2 K o = a(m d + m s )m u , (20) 
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where a = 3A/7T 2 /! 2 . It is only the baryon number Goldstone boson that remains 
massless. This mode is responsible for the superfluid nature of this phase (|120j) . 

Typically the pairing contribution to the pressure is small, being ~ AV be- 
cause only particles in a shell of width A at the Fermi surface pair. We can 
supplement the bag model pressure with the pairing contribution due to super- 
conductivity to obtain the pressure of the CFL phase (|122|) 

^CFL = ^ ~ ^2 + ~2 A V ~ B ■ ( 21 ) 

From this equation we can infer that, although the pairing contribution is para- 
metrically small by the factor A 2 //i 2 compared to the kinetic pressure, it could 
become the dominant contribution in the region where kinetic pressure and the 
bag pressure nearly cancel each other. Furthermore, because the quark-hadron 
phase transition is expected to occur in this regime, pairing will likely lower the 
quark-hadron phase transition density (|124II123[) . Similar trends are seen in other 
model calculations such as the Nambu-Jona-Lasino model of the quark EoS (125). 

Using Equation with a c g ~ 0.6 to describe the quark EoS and the APR 
model to describe the nuclear EoS, Alford et al. (126| find that the quark EoS 
and the nuclear EoS yield nearly the same pressure for a range of energy densities 
in the vicinity of nuclear density. This surprising finding, albeit within a specific 
model, indicates that a phase transition at these densities is unlikely to result 
in significant softening. Nonetheless it would have dramatic consequences for 
the transport and cooling properties because the excitation spectrum of the CFL 
phase is very distinct from that of nuclear matter. 

4.3.2 Role of the Strange Quark Mass It is a reasonable approxi- 
mation to neglect the up and down quark masses at densities of relevance to 
neutron stars where fi ~ 400 MeV. The strange quark mass m s ~ 200 MeV, how- 
ever, cannot be neglected. When the strange quark mass is large, comparable 
to the baryon chemical potential, the Fermi surfaces of the strange quarks and 
light quarks will be different. If this difference is similar to A, pairing involving 
strange quarks is suppressed. In the limit of infinite strange quark mass, i. e., 
in the absence of strange quarks, only the light quarks pair. This phase two- 
flavor superconducting phase is termed 2SC and is also characterized by pairs 
that are antisymmetric in flavor (118). Antisymmetry in color space excludes 
one of the three colors from participating in the condensation. Thus SU (3) co ior 
is broken down to SU (2) co i or and quarks of a particular color and three gluons 
remain massless. Furthermore, unlike in the CFL phase no global symmetries are 
broken. 

Early attempts to bridge the gap between the 2SC and the CFL phases can 
be found in References (|127|) and ((128). The authors found that the CFL pairing 
scheme was preserved when m s ~ y/2/iA, whereas for larger m s , a first-order 
transition to the 2SC phase occurs. Bedaque & Schafer showed that the strange 
quark mass appeared in the effective theory for Goldstone bosons in the form 
of a chemical potential for antistrangeness (|129j) . When the chemical potential, 
fis ~ m 2 /2/i, exceeds the mass of K° which is lightest meson with antistrangeness, 
these bosons condense in the ground state. In the CFL-fT phase hypercharge or 
strangeness symmetry is spontaneously broken resulting in the appearance of 
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a massless Goldstone boson. For a detailed discussion of these novel meson- 
condensed phases and their role in the birth and evolution of neutron stars see 
Reference (fT30"|) . 

A finite m s which is larger than the up and down quark masses leads to a 
suppression in their number and requires electrons to ensure neutrality. The 
electron chemical potential needed to accomplish this is fi e ~ m^/4/i. This in 
turn results in a splitting between the up and down Fermi levels because in beta- 
equilibrium n u = fj, — 2fi e /3 and (jl^ = \i + /x e /3. In two-flavor quark matter 
charge neutrality requires a relatively large \i e ~ \x because strange quarks are 
either absent or too few in number to furnish the necessary negative charge. 
Consequently, pairing will be favored only in strong coupling (A ~ \i) and in 
general this disfavors the 2SC phase in the neutron star context (|131|) . However, 
in Nambu-Jona-Lasino model calculations of quark matter one typically finds 
that the effective strange quark mass is large and the quark pairing strength is 
strong. Both of these effects favor the possibility of realizing the 2SC phase at 
intermediate densities of relevance to compact stars (|132[) . If, however, A <C ^ 
and pairing between up and down is not possible, then quarks of the same flavor 
can pair in the spin-one channel. The spin-one phase is characterized by a locking 
between spin and color, and model predictions indicate a small gap of the order 
of one MeV 

BCS pairing between two species of fermions results in a locking of the Fermi 
surfaces in momentum space and hence enforces equality in the numbers. This 
rigidity arises owing to the finite gap which requires a finite energy to produce a 
quasi-particle excitation. For example, in the CFL phase this locking ensures 
equal number densities for all three quark flavors even though m s > m^ght 
(|134U135f) . In the 2SC phase, the number densities of the up and down quarks are 
equal even when \i e ^ 0. This persists until fx e ~ 2A when it becomes possible 
to generate a number asymmetry at the expense of pairing. However, in simple 
systems, without long-range interactions, as fj, e increases one encounters a first- 
order phase transition from the BCS phase to the normal phase when /i e > \/2A 
. Surprisingly, even in weak coupling the phase structure of asymmetric fermion 
systems is not well known. Several candidate ground states have been suggested. 
They include(a) the mixed phase, where normal and superconducting phases co- 
exist (|136f) (b) the LOFF phase (named after the authors Larkin, Ovchinnikov, 
Fulde, and Ferrell who first suggested this possibility (|137|) ). where the gap ac- 
quires a spatial variation on a length scale ~ 1/A and (c) the gapless superfluid 
phase (also known as the breached-pair phase) which is a homogeneous super- 
fluid phase with gapless fermionic excitations due to a finite occupation of quasi- 
particle levels H138|) . Recent work has shown that in systems with short-range 
forces the gapless superfluid phase is unstable with respect to phase separation 
(136). However, as we discuss below, all of these possibilities could be relevant 
in dense quark matter. 

The presence of long-range forces and resulting constraint of charge neutrality 
in the bulk influences the competition between the various proposed phases of 
asymmetric Fermi systems. For example although the mixed phase is favored 
over the gapless phase in neutral systems the neutrality constraint can stabilize 
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the gapless phase with respect to phase separation in two-flavor (|139j) and three- 
flavor (|140|) quark matter. These phases which are termed the gapless 2SC and 
gapless CFL phases are remarkable because they are superconducting and yet 
have a large number of gapless excitations near the Fermi surface. In strong 
coupling, quantum Monte Carlo results show that gapless phases are stable with 
respect to phase separation Q141|) . However, in weak coupling there remain several 
questions that need to be addressed to ascertain if these phases or minor variants 
of them are indeed stable. In particular in QCD the magnetic gluon masses 
are imaginary in the weak coupling gapless phases. This signals an instability 
known as the chromomagnetic instability and occurs both in gapless 2SC (|142|) 
and gapless CFL phase (|143|) . Because the Meissner masses can be related to 
the gradient energy associated with spatial deformation of the pairing field some 
type of heterogeneity is likely and it appears that either the mixed phase or the 
LOFF phase is favored over gapless 2SC (|144j) . This may also be the case for 
three-flavor quark matter Q145|) . 

4.3.3 Absolutely Stable Strange Quark Matter The conjecture 
that matter containing strange quarks could be absolutely stable is several decades 
old fl_46J . This could occur if SQM has a lower energy per baryon at zero pressure 
relative to heavy nuclei such as 56 Fe, for which the energy per baryon is ~ 930 
MeV. If this were true, nuclear matter and nuclei would only be metastable, al- 
beit with a very large lifetime because it would take several simultaneous weak 
interactions to generate the strangeness needed to lower the energy. However, 
this situation may not be difficult to achieve in the hot and dense interior of a 
neutron star. Once strangeness is seeded somewhere in the star it should be able 
to catalyze the conversion of the entire star on a short time scale, producing a 
"strange star", i.e., a compact star made entirely of SQM. For a recent review of 
SQM and its role in compact stars see Reference (|147|) . 

We briefly mention a few recent findings of interest. For SQM to be preferred 
over nuclei or nuclear matter its pressure must be larger than the pressure in 
nuclear matter when the quark chemical potential fi ~ 930/3. From Eauation l21l 
we see that in the Bag model this would require a relatively small bag constant. 
However, the bag constant cannot be too small otherwise even two-flavor quark 
matter would have alower energy than nuclei and these would instantly convert to 
quark matter as it would not require the simultaneous weak interactions needed 
to produce strange quarks. In a simple model where a e fj ~ 1, A = and m s = 150 
MeV a bag constant given by B 1 / 4 ~ 145 MeV predicts SQM to be absolutely 
stable while ensuring that the 2 flavor quark matter is unstable. A large gap of 
the order of 100 MeV has the effect of increasing the effective bag constant at 
which quark matter can become stable. Recent work on the effects pairing in 
SQM can be found in Ref. (tT21|l . 

Early work indicated that SQM would be homogeneous liquid at arbitrarily low 
pressure. This is unlike low-density nucleon matter made of nuclei and electrons 
which is a heterogeneous solid with nuclei embedded in a background electron 
gas. Reference (|148j) investigated the possibility that SQM could also be a het- 
erogeneous solid at low pressure. Using a model-independent analysis this work 
shows that at relatively low density SQM could also become heterogeneous and 
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phase-separate into positively charged nuggets of quark matter embedded in a 
background electron gas - much like matter with nuclei at low pressure. With 
increasing pressure the quark phase would be populated with voids of electrons. 
If stable, the heterogeneous crust could have important consequences for surfaces 
of strange quark stars. Unlike the conventional scenario for strange stars which 
are characterized by an enormous density gradient and a large electric field (|149[) 
the heterogeneous solid crust would occupy a large radial extent (Ai? ~ 50 — 100 
m) and obviate the need for an electric field. 

4.3.4 Excitation spectrum and response properties Response prop- 
erties such as neutrino emissivities, viscosity and the low-temperature specific 
heat have been calculated in both normal and superconducting phases of quark 
matter. In the normal phase, a large number of particle-hole excitations exist 
at the Fermi surfaces of the quarks and result in enhanced response and a large 
specific heat. Consequently neutrino cooling due to processes analogous to the 
direct URCA reaction but involving quarks is rapid ({150p. The bulk viscosity 
is also enhanced owing to the weak interaction rate that converts d «-> s quarks 
(|151[) . In contrast, in the CFL phase where all quarks are gapped, the low-lying 
excitation spectrum is sparse and only the baryon number Goldstone mode con- 
tributes at temperatures of relevance to neutron star cooling. Consequently the 
neutrino emissivity (152), specific heat and bulk viscosity |153J are all relatively 
negligible. This suppression in the neutrino processes also persists in the CFLi^ 
phase(;15|J). In less-symmetric superconducting phases such as 2SC where un- 
gapped quark excitations exist the situation is similar to that in the normal 
quark matter. 

The existence of either gapless color superconductivity or the heterogeneous 
phases such as the LOFF phase or the mixed phase can lead to potentially ob- 
servable consequences. In the gapless phase the number of gapless excitations is 
anomalously large. The resulting specific heat is a factor a/A/T larger than even 
the normal phase and could potentially impact neutron star cooling at late time 
(|155|) . The coexistence of heterogeneity and superfluidity is key to understand- 
ing glitches. In this sense the heterogeneous superconducting phases could play a 
role. These connections, between the properties of the phase and the observables, 
are still in their infancy and warrant further work. 

4.4 Nature of the Phase Transition 

The novel phases discussed above could occur either via a first- or second-order 
phase transition. In first-order transitions phase coexistence is possible. A het- 
erogeneous phase with coexisting phases is termed the mixed phase. Here charge 
neutrality is enforced globally. Each of the two coexisting phase have opposite 
electric charge and the volume fractions adjust to ensure overall neutrality. 

Heterogeneous phases are commonplace in the terrestrial context. For example 
a solid can be viewed as a mixed phase composed of positively nuclear matter 
(residing inside nuclei) coexisting with a negatively charged electron gas phase. In 
the relatively low-density region of the neutron star crust we encounter a similar 
mixed phase that upon increasing density changes to accommodate neutrons in 
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Figure 13: Behavior of the electron chemical potential and the pressure of homo- 
geneous neutral nuclear and CFL matter and of the mixed phase, all as a function 
of the quark chemical potential \i. 

the electron gas phase. The physics of the neutron-rich mixed phase in the crust 
is reviewed in detail in Ref. (|156j) . Glendenning noted that similar considerations 
apply to first-order transitions in the context of high-density phase transitions. 
This applies in general to all first-order phase transitions with two conserved 
charges ()157j) . In the neutron star context the conserved quantities are baryon 
number and electric charge. Unlike simple first-order transitions such the water- 
vapor transition which is characterized by one conserved charge (number of water 
molecules) where phase coexistence occurs at specific value of the pressure, the 
phase coexistence in dense matter occurs over a finite interval in pressure owing 
to the presence of the extra degree of freedom, namely electric charge. 

To illustrate the physics of first-order phase transitions and the role of surface 
and Coulomb energies in the mixed phase we consider an explicit example. The 
phase transition from nuclear matter to CFL quark matter is a first-order phase 
transition. The nuclear phase has no strangeness and the bulk quark phase 
has no electrons. The possibility of phase coexistence between these phases was 
investigated in Ref. (122). We highlight some of the main findings here. In 
Figure 1131 the pressure of the bulk nuclear, bulk CFL and mixed phases are 
shown as a function of /j,, the quark chemical potential. At intermediate values 
of /i, the mixed phase has larger pressure and is therefore favored over both the 
nuclear and CFL bulk phases. The electron chemical potential, /i e , required to 
ensure charge neutrality in the bulk nuclear phase, which grows with ^l, is shown. 
In the mixed phase neutrality requires a positively charged nuclear phase and a 
negatively charged CFL phase. This is easily accomplished by lowering /i e . The 
decreasing /U e in the mixed phase shown in the figure is obtained by requiring 
equal pressures in both phases at a given /i. 
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In the mixed phase the Coulomb and surface energy costs must be met. The 
results shown in Figure El ignored these corrections. Only bulk free energy is 
included; surface and Coulomb energy is neglected. The mixed phase occurs 
between A and D. The vertical line connecting B and C denotes the /x at which 
the pressures of neutral CFL and nuclear matter are equal. This is where a sharp 
interface may occur. The pressure of the mixed phase exceeds that of neutral 
CFL or neutral nuclear matter between A and D. Were this the whole story, the 
mixed phase would evidently be favored over the sharp interface. 

In a simple description of the mixed phase one considers a thin boundary 
between the coexisting phases. A unit cell of the mixed phase is defined as the 
minimum size region that is electrically neutral. Three different simple geometries 
are considered: spheres (d = 3), rods {d = 2) and slabs {d = 1) (158). In each of 
these cases, the surface and Coulomb energy cost are given by 

E s = d X aQCD , (22) 

E c = 2vr a cm f d (x) (AQ) 2 r 2 , (23) 

where d is the dimensionality of the structure (d = 1,2, and 3 correspond to 
Wigner-Seitz cells describing slab, rod and droplet configurations, respectively), 
a is the surface tension and AQ = (Jnuciear — QcFL+kaons is the charge-density 
contrast between the two phases. The other factors appearing in Equations (|22jl 
and (|23|) are x, the fraction of the rarer phase that is equal to x where \ ^ 0.5 
and 1 — x where 0.5 ^ % ^ 1; ro, the radius of the rarer phase (radius of drops or 
rods and half-thickness of slabs); and f d (x), the geometrical factor of order one 
that arises in the calculation of the Coulomb energy (156). The first step in the 
calculation is to evaluate ro by minimizing the sum of E and E s . The result is 

7 -i 1/3 

d x a QCD 1 (24) 



4vr a em f d (x) (AQ) 2 



Using this value of ro in Equations (|22|) . (|23j) the surface and Coulomb energy 
cost per unit volume is obtained 

ES + EP = 1 (4tt a cm d 2 f d (x) x 2 ) 1/3 (AQ) 2 / 3 a 2 £ D . (25) 

The lowest curve in Figure ITU shows Afi, the difference between the free-energy 
density of the mixed phase (calculated without including the surface and Coulomb 
energy cost) and the homogeneous electrically neutral nuclear and CFL phases 
separated by a single sharp interface, whichever of the two is lower. The mixed 
phase has lower bulk free energy, so AJ7, plotted in Figure ITU is negative. The 
remaining curves in Figure El show the sum of the bulk free-energy difference 
A$7 and (E s + E C ), and the surface and Coulomb energy cost of the mixed phase 
calculated using Equation El for droplets, rods and slabs for three different values 
of (Tqcd- Careful inspection of the figure reveals that for any value of ctqcd, the 
mixed phase is described as a function of increasing density by a progression 
from drops to rods to slabs of CFL matter within nuclear matter to slabs to rods 
to drops of nuclear matter within CFL matter. This is the same progression of 
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Figure 14: The free energy difference between the mixed phase and the homo- 
geneous neutral nuclear and CFL phases. In the lowest curve, the surface and 
Coulomb energy costs of the mixed phase are neglected, and the mixed phase 
therefore has the lower free energy. Other curves include surface and Coulomb 
energy for different values of oqcd and different mixed phase geometry. 



geometries seen in the inner crust of a neutron star (|158f) or in the mixed phase at 
a first-order phase transition between nuclear matter and unpaired quark matter 
(j!59j) or hadronic kaon condensate (|160|) . We have also checked that for oqcd = 
10 and 40 MeV/fm 2 , with the mixed-phase geometry at any x taken to be favored, 
the regions of both the rarer and more common phases (ro and its suitably defined 
counterpart) are always less than 5 — 6 fm. In general, uniform regions of charge 
can exist only on length scales small compared to the Debye screen length. The 
Debye screening length in the quark and hadronic phase are typically 5 — 10 
fm. When the size of the charged regions becomes comparable to the Debye 
screening length it becomes important to account for spatial variations of the 
charge density. This will influence the surface and Coulomb energy estimates 
presented in Equation 1251 A detailed discussion of the importance of these finite 
size effects is presented in References (|161|) and (jl62 ). 

For any given oqcd 5 the mixed phase has lower free energy than homoge- 
neous neutral CFL or nuclear matter whenever one of the curves in Figure El 
for that oqcd is negative. We see that much of the mixed phase will survive 
if ctqcd — 10 MeV/fm 2 , whereas for oqcd ^ 40 MeV/fm 2 the mixed phase 
is not favored for any \x. This means that if the QCD-scale surface tension 
ctqcd ^ 40 MeV/fm 2 , a single sharp interface will be favored. The interface 
is characterized by a bipolar charge distribution, resulting in an intense electric 
field that ensures that the electric charge chemical potential is continuous across 
it (see Ref. (|T22l for details). 
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5 CONFRONTING OBSERVABLES 

Observations relating to the structure and the evolution of compact stars fur- 
nish complementary information about the interior state. The former provides 
information about the EoS of dense matter and probes the high-energy or short- 
distance aspects, whereas the latter provides information about the low- lying 
excitation spectrum. This complementarity potentially allows us to probe the 
phase structure of matter at extreme density, as phase transitions can drastically 
alter the low-lying spectrum without strongly influencing the EoS or vice versa. 
For example, nucleon superfluidity and superconductivity dramatically alter the 
excitation spectrum (exponentially suppressing neutrino emission rates) but play 
no role in the EoS. In contrast, hyperons strongly soften the EoS but may not 
affect the linear response properties to the same degree. 

5.1 Bulk Properties and Structure 

Mass and radius are certainly the most natural observables that probe the EoS of 
the dense interior. Given the relationship between pressure and energy density, 
the general relativistic equation of hydrostatic equilibrium (termed the Tolman- 
Oppenheimer-Volkoff, TOV, equation) uniquely determines the structure of the 
star, in particular its mass and radius, for a given central density (85). For 
nonrotating stars, each high-density EoS uniquely specifies a mass-radius curve. 
In the context of neutron star structure, EoSs are characterized as soft or stiff 
relative to each other on the basis of the ratio of pressure to energy density. This, 
however, depends on the density at which the comparison is made. The EoS with 
the larger (smaller) ratio at a specific density is termed stiff (soft). An EoS that 
is soft on average will lead to more compact stars, compared with an EoS that 
is stiff, and reach a smaller maximum mass M max . This is easily understood by 
noting that the energy density is the source of gravity, whereas pressure provides 
resistance to the gravitational squeeze. 

From model studies it is empirically known that the maximum allowed mass 
is most sensitive to the EoS at the highest (supranuclear) density, whereas the 
radius is sensitive to the EoS in the vicinity of nuclear density fl63f). In particular, 
for stars composed of nucleon matter, the radius appears to be fairly sensitive to 
the density dependence of the nuclear symmetry energy. The symmetry energy 
can potentially be inferred from terrestrial measurements of the neutron-skin 
thickness (164J and probably from heavy- ion experiments (|165[) . For a recent 
review of the role of the nuclear symmetry energy in neutron star structure and 
terrestrial experiment, see Reference (166:). These experiments likely constrain 
models of the nuclear EoS that currently differ significantly at both low and high 
density. For example, the mean-field models 1)85(1 tend to be stiff at low density 
and soft at high density compared with variational calculations such as those 
obtained by Akmal et al. (J72J). 

Phase transitions to novel phases at supranuclear density typically result in a 
softening of the EoS. This arises because these transitions furnish new degrees of 
freedom that contribute more to the energy than to the pressure. For example, a A 
hyperon that replaces a neutron at the Fermi surface has a much lower momentum 
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Figure 15: Mass-radius diagram: The three generic classes of equations of state 
are depicted with observational constraints. 

and a larger rest mass. In general, this is true for all hadronic transitions studied 
to date. Clearly the extent of the softening depends on the details of both the 
nuclear and exotic EoSs. To successfully employ a mass-radius constraint that 
infers the phase structure of the dense interior, softening should be significant. 
This quantitative expectation is, however, difficult to assert and may not be 
true in general. For example, early studies that employed the naive bag model 
indicated that the quark-hadron transition would lead to significant softening, but 
recent work, which takes into account corrections due to quark-quark correlations 
and superconductivity, suggests that the quark EoS might mimic nuclear behavior 

dim). 

Figure Hoi depicts the M — R relationship of three generic classes of stars: (a) 
nuclear stars, (b) hybrid stars, and (c) strange (SQM) stars. Here we show bands 
rather than individual M — R curves for specific models to indicate the inherent 
model dependence of these theoretical predictions. As discussed above, the broad 
range of predicted radii for nucleon EoSs will be narrowed in the near future 
owing to neutron-skin thickness and probably also to heavy- ion experiments. As 
anticipated in hybrid stars, where there is a phase transition from a nuclear to an 
exotic phase at supranuclear density, softening leads to smaller radii and lower 
maximum masses. In contrast, the most exotic scenario involving absolutely 
stable SQM overlaps the regions accessed by nuclear and hybrid stars for masses 
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in the observed range (1-2 M & ). Numerous constraints on the EoS are depicted 
on the figure. General relativity requires that the radius be larger than the 
Schwarzschild radius Rs and causality of the EoS requires that the speed of 
sound c s = {dP/dp) 1 / 2 be smaller than the speed of light. This latter constraint 
is approximately equivalent to i? max > l-5-Rsmax> where i? max and Rs max are 
the radius and the Schwarzschild radius, respectively, of the maximum mass 
star {167 ). (Note that the softest nuclear EoSs, resulting from nonrelativistic 
models, do violate causality at densities close to the maximum density reached 
at M max .) The other curves in Figure El refer to various constraints arising from 
observations: (a) lower limits on the mass of PSR J0751+1807, at the 68% and 
90% C.L. (3), and lower limits on the estimate (68% C.L.) of the mass of 4U 1700- 
37 (|14|) . which may be a neutron star or a black hole; (b) the mass-shedding limit 
on the maximum radius of the fastest known pulsar PSR J1748-2446ad (|69|): (c) 
90% C.L. of the radius at infinity Roc of the quiescent LMXB X7 (|44j) discussed in 
Section ^. 2, 41 and finally (d) the minimal radius of glitching pulsars deduced from 
the constraint I s /I c > 1.4%, (see Section 12.4(1 . assuming the glitch reservoir is 
the neutron superfluid in the crust. Converting this moment of inertia constraint 
into a mass-radius constraint simply requires integrating the TOV equation from 
the crust-core boundary to the surface: The pressure Pt at this transition point 
is, however, uncertain, with values ranging from 0.25 up to 0.65 MeV fur 3 (68) 
and the two I s /I c curve in Figure El labelled 0.25 and 0.65 correspond to these 
two extremes. For each of the two classes of EoSs, nucleon or hybrid, low P t 
EoSs give small radii at low masses, whereas high Pt give large radii: The glitch 
constraint is hence much more restrictive for EoSs which are soft at densities 
around nuclear matter density. 

These various constraints shown in Figure El are already sufficient to draw 
some conclusions about the high-density EoS: 

1. The softest range of the hybrid EoSs shown in the figure has a maximum 
mass above the accuretely measured mass of the PSR B1913+16 compan- 
ion, 1.4408 ± O.OOO3M (fTfiKjt . but they are now ruled out at the 90% C.L. 
by the new measurement of PSR J0751+1807. Moreover, all hybrid EoSs 
are ruled out at the 68% C.L. by the same measurement. 

2. Hybrid EoSs which are soft at low density, i.e., have a low Pt ~ 0.25 MeV 
fm~ 3 , are also strongly disfavored by the glitch I s /I c ~ 1-4% constraint. 

3. The same EoSs are also in strong disagreement with the Roo measurement 
of X7. More reliable measurements of neutron star masses in quiescent 
LMXBs are needed to confirm this conclusion and are likely to become 
available in the near future, as we discussed in Section 12.2.41 

4. The softest nucleon EoSs, which have Pt ~ 20-30 MeV fm -3 , are marginally 
consistent with the glitch constraint and need glitching pulsars to have 
masses M ~ 1. l.U. . 
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Figure 16: Slow, "standard", cooling via the modified URCA process versus 
fast neutrino cooling from three neutrino emission processes with emissivities 
e u = 10 n x Tg erg cm -3 s~ 1 . Trajectories marked as N assume neutrons in the 
core are normal, whereas cases marked SF assume neutrons are paired in the 
3 P2 channel in the whole core with a T c of the order of 2-3 x 10 9 K. The fast- 
cooling case with n=24, 25, and 26 corresponds approximately to the emissivities 
of a kaon condensate, a pion condensate, and the nucleon direct URCA process, 
respectively. Note that at early times, t ~ 30 yrs, as long as the age of the star 
is inferior to the thermal diffusion time scale through the crust, all models result 
in the same surface temperature and luminosity. Neutrino emission from the 
Cooper pair formation process is not taken into account in these simple models. 
Figure adapted from Reference (170). 



5. Given the broad range of predicted radii for strange stars, none of the 
present constraints can exclude them or make a strong case for their exis- 
tence. 

5.2 Low-Lying Spectrum and Thermal Behavior 

In contrast to bulk properties, the response properties influence observable as- 
pects such as thermal evolution and typically probe energy scales of the order 
of tens to hundreds of keV, and are very sensitive to the low-lying excitation 
spectrum. In some respects, they also tell us about the subtle properties that 
determine the phase structure of dense matter. 

5.2.1 The Cooling of Isolated Neutron stars The long-term cooling 
of neutron stars is the best studied example in this class. This topic was recently 
reviewed in detail (|169U31[) and we present only a short description here. This 
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cooling is driven by neutrino emission from the core of the star when its tem- 
perature is above 10 7 — 10 8 K or 1 — 10 keV. The time scale for cooling can be 
estimated in terms of two key microscopic ingredients, namely the neutrino lumi- 
nosity, L u , and the total specific heat, C v , supplemented by the surface photon 
luminosity, L 7 , from a simple energy balance consideration as follows 

^L-C — - -L -L (26) 
dt ~ Lv dt ~ Lu ^ ' {2b) 

where E t h is the thermal energy content of the star. For a degenerate Fermi 
system, the specific heat C v is linear in T, i.e., C v ~ C -T, where the constant C 
depends on the details of the star's structure. The neutrino luminosities can be 
split, in a first approximation, into fast or slow processes and can be written as 
L fast ~ 7v f • T 6 or L" low ~ iV s • T 8 , respectively. Finally, L 7 oc T e 4 and because, 
very roughly, T e oc T 1 / 2 , one can also write L 7 ~ S ■ T 2 . Owing to their very 
different T dependences, at early times L v 3> L~, whereas when T is sufficiently 
low the converse is true. From Equation 1261 the two neutrino cooling time scales 
for fast and slow neutrino emission are 

fast C _ 4 minutes 
4iV f T 4 ~ Tg 

slow _ C „ 6 months 
6N s T 6 ~ T 6 > 

respectively. In obtaining the approximate numerical estimates (4 minutes and 
6 months) we have used typical values of the constants C, N , and iV s , which 
are characteristic of degenerate dense hadronic matter (|31|). For fast cooling, 
L u dominates L 7 until T ~ 10 6 K, whereas for slow cooling this happens when 
T ~ 10 8 K. Coincidentally, in both cases this turns out to occur at ages of the 
order of 10 5 years. 

In general, exotic phases almost invariably result in fast neutrino emission (the 
CFL phase is the only known exception), and detection of a fast-cooling neutron 
star would be an argument in favor of the presence of one of these phases. How- 
ever, because nuclear matter with a high enough proton fraction also results in 
fast neutrino emission through the nucleon direct URCA process, such a detec- 
tion would only be an indication and not proof of one of these phases. Moreover, 
because pairing results in a strong suppression of the neutrino emissivity of all 
processes in which the paired component participates, an exotic phase may be 
present and may not manifest in the thermal evolution of the star. We illustrate 
these considerations in Figure 1161 in which the time evolution of the thermal lu- 
minosity Lqo °f a model neutron star for the case of the slow "standard" cooling 
from the modified URCA process is compared with three different fast neutrino 
emission processes. Each scenario is presented in two forms, one in which neu- 
trons are normal and one in which where they are assumed to be superfluid in 
the whole core. The occurrence of fast cooling leads to thermal luminosities, at 
ages ~ 10 2 -10 5 yrs, which are about three orders of magnitude lower than in 
"standard" cooling. However, when the suppression of the neutrino emissivity 
by pairing is taken into account, the differences can be much smaller. Higher 
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Figure 17: Comparison of the prediction of the minimal cooling paradigm with the 
observational data presented in Figures [7] and |HJ Figure adapted from Reference 

(El. 



values of T c can even render fast and slow cooling indistinguishable (|171j) . More- 
over, once the fast neutrino emission is controlled by pairing gaps, the difference 
between cooling efficiencies of various fast processes becomes much smaller than 
when they proceed unhindered, as the three fast-cooling superfluid (SF) curves 
of Figure ED show when compared with corresponding curves with normal (N) 
neutrons. Considering the plethora of channels for fast neutrino emission when 
the number of degree of freedom increases, (a) pion or kaon condensates, (b) hy- 
perons, and (c) deconfined quark matter, it appears unlikely that these kinds of 
studies by themselves will be able to distinguish between types of exotica. Ref- 
erence (|172j) explicitely showed that model stars with only nucleons, or nucleons 
with hyperons, or nucleons with deconfined quark matter or nucleons with hyper- 
ons and deconfined quark matter can lead to essentially indistinguishable cooling 
trajectories once uncertainties about the size of the various pairing gaps are taken 
into account. With respect to strange stars, if the SQM surface is covered by a 
thin envelope of normal nuclear matter, their thermal evolution would be similar 
to the one of the more mundane neutron stars (|31U173f) and only in the case of a 
bare quark surface would their thermal evolution be drastically different (|174|) . 

Given these intricacies, Page et al. (33) introduced the minimal cooling sce- 
nario as a paradigm in which all exotic phases are assumed to be absent and the 
proton fraction is also assumed to be below the direct URCA threshold. This 
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paradigm is a benchmark against which data are compared so that observations 
that cannot be accommodated within the predictions of this scenario are serious 
evidence for the presence of some fast-cooling agent. Minimal cooling is not naive 
cooling: It considers all others processes that can affect the thermal evolution of 
the star and are based on very standard physics. Two standard ingredients can 
influence on the star's cooling significantly. The first is the occurrence of neutron 
superfluidity and/or proton superconductivity, which suppresses the specific heat 
and the neutrino emission by the modified URCA and bremsstrahlung processes 
but also opens the new neutrino emission channels owing to the constant forma- 
tion and breaking of Cooper pairs. The net effect of pairing can hence be either 
a reduction or an enhancement of Lf ow . The second important ingredient is the 
chemical composition of the uppermost layers of the star, the envelope, which 
controls the photon luminosity L 7 : A light-element envelope, H, He, C, or O, is 
less insulating than a heavy-element, iron- like envelope and results in a larger L 7 
and T e for a given internal temperature. A light-element envelope has the effect 
of making the star brighter during the neutrino cooling era because at that time 
L 7 is negligible compared to L u , and the surface temperature T e just follows the 
evolution of the internal temperature while accelerating the cooling owing to its 
higher L 7 during the photon cooling era. In Figure 1171 the predictions of the 
minimal cooling paradigm are displayed and compared with the observational 
data presented in Section [2.2.31 Models with a light-element and heavy-element 
envelope are separated, and for each class the spread in prediction is owing to 
the uncertainty in the size of the neutron and proton pairing gaps. Note that the 
luminosities at ages 10 2 -10 5 yrs can be quite lower than those of the "standard" 
cooling of Figure ED because of the enhanced neutrino emission from the Cooper 
pair formation process. Overall the agreement between the theory and the ob- 
servational data is quite good. Two objects, PSR J0205+6449 (in SNR 3C58) 
and RX J0007.0+7302, nevertheless have upper limits below these predictions 
and are good candidates for the occurrence of fast neutrino emission, controlled 
by pairing, but they may also be interpreted as evidence for medium-enhanced 
modified URCA process (176). The four upper limits, marked a, b, c, and d in 
FigureEl from Kaplan et al.'s (|3*5|l search, would, however, be definitive evidence 
for fast cooling if they can be proven to refer to neutron stars and not black holes. 
Similar conclusions have been reached by the St. Petersburg group (|175[) . which 
has developed a complementary version of minimal cooling. 

5.2.2 The Thermal Behavior of Neutron Stars in Low-Mass X-Ray 
Binaries We presented in Section \2 .31 observations of neutron stars in LMXBs 
that are undergoing transient phases of accretion separated by periods of quies- 
cence. Brown et al.'s (|49|) mechanism of heating of the star provides a simple 
explanation for the hot thermal spectrum observed in quiescence in the cases 
of Aql X-l, 4U 1608-522, and CX 1, but the results of Table □ clearly show 
that many systems have quiescent luminosities much below the prediction of 
this model. A natural explanation for this discrepancy is the presence of fast 
neutrino cooling in the core of these too-cold neutron stars (62). This result 
condradicts sharply with the fact that most isolated cooling neutron stars have 
thermal luminosities in agreement with the results of the minimal cooling see- 
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nario, i.e., they show no conclusive evidence for the occurrence of fast neutrino 
emission. However, remembering that LMXBs are very long-lived systems in 
which the neutron star may accrete a significant part of a solar mass, these cold 
neutron stars may simply be massive enough to be above the threshold for fast 
neutrino emission, whereas most of the isolated cooling neutron stars are below 
this threshold. These results would be in good agreement with the expected mass 
range of neutron stars emerging from core collapse supernovae as illustrated in 
Figure 01 and are corroborated by the measured masses of PSRs in PSR+NS sys- 
tems exhibited in Figure [2 whereas the measured masses of PSR+WD systems, 
offsprings of LMXBs, confirm the existence of more massive neutron stars. 

5.3 Emerging Trends and Future Prospects 

5.3.1 Limiting Spin Frequency The first, absolute, limiting frequency 
for a fast-spinning pulsar is the mass-shedding limit Oms > which we used as 
a constraint in Figure 1151 from the fastest known pulsar PSR J1748-2446ad. 
However, gravitational radiation reaction-driven instabilities such as the r-mode 
instability can result in a more stringent limit on the spin frequency. These 
modes generate gravity waves (which dissipate angular momentum) that, instead 
of damping, amplify mode amplitude (|177f) . Viscosity is then the main source 
of damping and, depending on the bulk and shear viscosity of dense matter, 
different limiting spin frequencies are possible. At very low temperature, damping 
is due to the shear viscosity, whereas at very high temperature, it is due to the 
bulk viscosity. For these extreme temperatures, r-modes are effectively damped. 
At intermediate values of temperature, depending on the bulk viscosity and its 
temperature dependence, a limiting frequency as low as 50% of f^MS is possible. In 
LMXBs where the neutron stars are spun up to millisecond periods, a population 
study shows a sharp drop in the pulsars with spin frequency O ~ 730 Hz ([T7R|) . If 
this limiting behavior is indeed due to gravity waves from the r-mode instability, 
it provides useful information about the viscous damping rate in the star, which 
can be calculated for different EoSs for the neutron star interior. For example, 
a strange quark star made entirely of CFL quark matter would be incompatible 
with these observations because all low-energy weak interaction processes that can 
contribute to the bulk viscosity are exponentially suppressed by exp(— A/Zc^T), 
where A is the pairing gap and T is the temperature (179). 

5.3.2 Galactic Supernova Neutrinos The enormous gravitational bind- 
ing energy, Eq ~ GM 2 /R ~ 3x 10 53 ergs, gained during a core collapse supernova 
is stored inside the newly born proto-neutron star as thermal energy of the mat- 
ter components and thermal and degeneracy energy of neutrinos. This energy 
then leaks out on a time scale determined by the rate of diffusion of neutri- 
nos in the dense core (180). The detection of ~ 20 neutrinos from SN1987A 
is testimony to this theoretical expectation. Current estimates indicate that we 
should see ~ 10, 000 events in SuperKamiokande and ~ 1000 events in the Sud- 
bury Neutrino Observatory (SNO) from a supernova at the center of our galaxy 
(distance=8.5 kpc) (|181[) . Detection would provide detailed information about 
the temporal structure of the neutrino signal and thereby probe the properties of 
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the dense inner core. Preliminary studies of the role of phase transitions in the 
evolution of proto-neutron stars and the supernova neutrino signal exist and show 
that the late time signal is sensitive to the high-density physics. In particular, a 
dramatic consequence of the transition to a very soft high-density EoS would be 
a delayed collapse to a black hole on the neutrino diffusion time scale (|182[) . Less 
dramatic changes to the neutrino luminosity can be detected for a galactic super- 
nova and could potentially probe significant changes to the neutrino mean free 
path that occur as a result of a phase transition. Like neutrino emissivities, the 
neutrino opacity of dense matter is sensitive to the phase structure and low-lying 
excitation spectrum of dense matter (|183|) . 

5.3.3 Gravity Waves The detection of gravity waves (|184j) from a bi- 
nary system in Laser Interferometer Gravitational- Wave Observatory-like detec- 
tors, which can measure frequencies up to 1 kHz, can measure the chirp mass 
M ch = fi 3 / 5 Mj /5 , where is the reduced mass and Mt is the total mass of the 
binary (|185|) . Perhaps more interestingly, during the inspiral the compactness 
M/R is also accessible and will thereby provide an estimate of the star radius 
(186). In addition, a first-order phase transition occurring in the inner core of a 
neutron star could lead to a density discontinuity if the surface tension is large. 
In Reference (|187|) . the authors show that this would affect the frequency spec- 
trum of the non-radial oscillation modes in two ways. First, it would produce 
a softening of the EoS, leading to more compact equilibrium configurations and 
changing the frequency of the fundamental and pressure modes of the neutron 
star. Second, a new nonzero-frequency g-mode would appear in association with 
each discontinuity. If these modes are excited in either a phase transition-induced 
mini-collapse or a glitch (binary neutron star mergers are unlikely to excite these 
high-frequency modes), they may be observable in the next-generation gravity 
wave detectors. The role of the EoS in binary neutron star mergers is also po- 
tentially interesting and largely unexplored. Reference (188 ) finds that if there 
is time for stable mass transfer to occur during the merger, it could provide in- 
formation about the mass-radius diagram and potentially distinguish between 
normal stars and strange quark stars. 

6 CONCLUSIONS 

To infer the phase structure and properties of matter at extreme density through 
observations of compact stars is a compelling but difficult task. Despite the 
rather rudimentary state of quantitative theoretical calculations in describing 
dense matter, the observational data are already providing guidance to and hint- 
ing at possible tensions between different models. In particular the prospect for 
a firm mass measurement of a heavy neutron star (or a low-mass black hole) and 
accurate radius measurements are excellent. Such measurements will directly 
confirm or rule out the possibility for a strong phase transition to a soft EoS 
in the vicinity of nuclear density, and presently available values already seem to 
disfavor a soft EoS. Observations of cooling neutron stars and the realization that 
superfluidity and superconductivity play a dominant role in neutrino cooling of 
the core has sparked renewed interest in theoretical studies of pairing correlations 
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in dense matter. Furthermore, the connection between seemingly different ob- 
servables such as long-term neutron star cooling, quiescent luminosity observed 
in LMXBs, and superbursts are indicating new trends that point to a diversity 
in the thermal evolution of neutron stars, with growing evidence in favor of the 
occurence of some fast neutrino emitting process(es), at least in heavy, > 1.4M , 
neutron stars. These developments call for a more global approach in the the- 
ory and simulations to confront different observables within a unified framework 
so that the sensitivity to the underlying dense matter physics can be properly 
assessed. 

Theoretical progress in describing dense nuclear matter and new insights into 
using terrestrial experiments on nuclei to better constrain EoS models provide 
compelling reasons to expect a quantitative description of neutron-rich matter 
for the densities of interest. Throughout this review we highlighted important 
qualitative differences and possible similarities between nucleon matter and other 
exotic phases. Minimizing the theoretical uncertainty for the predictions of the 
nucleon matter is crucial to confront observables. For example, a precise deter- 
mination of the density dependence of the nuclear symmetry energy can reduce 
the uncertainty simultaneously in radii and cooling rates predicted for the dif- 
ferent nuclear EoSs. The situation with exotic states of matter is less clear. 
Although softening and enhanced cooling appear to be common trends, there are 
important exceptions and quantitative estimates are often very model dependent. 
Nonetheless, as mentioned above, radius measurement smaller than 10 km or the 
discovery of low-mass black hole candidate with mass ~ 2M would be very sug- 
gestive of a phase transition. Better experimental or theoretical constraints on 
the hyperon-nucleon and kaon-nucleon interactions are necessary to ascertain if 
novel hadronic phases are relevant. Quark matter remains an especially interest- 
ing possibility. Color superconductivity is relevant, especially for the calculation 
of neutrino processes and cooling. However, an understanding of how a finite 
strange quark mass affects the phase structure is needed because this can result 
in vastly different low energy properties. 
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